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abstracts 
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Chapter  I : 

Lipscomb  and  King's  hindered  rotor  expression 
was  applied  to  the  rigid  body  routine  of  M.J.  Bennett 
and  B.M.  Foxman.  The  programme  SFLS5  was  modified 
to  calculate  the  contribution  of  the  hindered  rotor 
to  the  structure  factors  and  to  do  a  full  matrix  least 
squares  refinement  of  the  input  parameters.  The 
programme  was  tested  on  several  sets  of  published 
X-ray  data.  The  results  showed  that  the  model  gives 
correct  bond  lengths  without  recourse  to  "riding 
corrections".  The  barrier  parameter  Bd  is  related 
to  the  RMS  angular  displacement  of  the  ring  but  is  not 
necessarily  a  measure  of  activation  energy. 

Chapter  II: 

The  compound  ferrocinium  picrate  [Fe  (C._H_ )  _  +  ] 

5  5  2 

[H^CgN^O^  ]  crystallizes  in  space  group  Cmcm  (Z  =  4) 

o  o 

with  cell  dimensions  a  =  12.513(7)A,  b  =  20.267(9)A, 

o 

c  =  6.903 (7) A.  The  structure  was  refined,  using  SFLS5 
as  modified  for  the  hindered  rotor  model,  to  a  conven¬ 
tional  R  factor  of  6.0%.  The  compound  consists  of  two 
ionic  species,  the  ferrocinium  cation,  and  the  picrate 
ion.  Both  ions  sit  on  sites  of  mm  symmetry  and  show 


some  disorder. 
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Chapter  III: 


Hydridotriphenylsilyl  (Tr-cyclopentadienyl) 

dicarbonyl  manganese,  (i-C  H  )  (OC)  0Mn  (H)  Si  (CrH_)  ,  , 

do  Z  6  5  3 

crystallizes  in  the  space  group  P2^/n  with  unit  cell 
dimensions  a  =  13.200(3)A,  b  =  17.328(4)A,  c  =  9.438  (2)A, 
3  =  92.47(4)°,  Z  =  4.  The  structure  was  refined  to  a 
conventional  R  factor  of  4.2%.  The  geometry  around  the 
manganese  atom  is  that  of  a  distorted  square  pyramid, 
with  that  of  the  silicon  atom  a  distorted  tetrahedron. 

The  hydridic  hydrogen  was  located  on  full  data  and  low 
angle  difference  maps  and  was  subject  to  refinement. 

It  is  located  in  a  bridging  position,  lying  1.54(4)A 

o 

from  the  manganese  and  1.76 (4) A  from  the  silicon. 

Chapter  IV: 

Dihydrooctacarbonyl  dirhenium,  H  ReJCO)., 

Z  Z  o 

o 

crystallizes  in  space  group  P21/c,  a  =  8.96 (2)A, 
b  =  11. 62 (2) A,  c  =  12 . 85 (2) A,  3  =  109.2(1)°,  with 
four  molecules  per  unit  cell.  The  structure  was 
refined  to  a  conventional  R  factor  of  10.5%,  all  atoms 
isotropic.  There  is  octahedral  geometry  about  each 
rhenium,  the  hydrogens  being  assumed  to  occupy  the 
vacant  co-ordinated  positions.  The  doubly  hydrogen 

o 

bridged  Re-Re  bond  of  2.896 (3) A  is  shorter  than  normal 
Re-Re  single  bonds. 


The  topics  covered  in  this  thesis  involve  X-ray 

diffraction.  Undefined  symbols  have  their  standard 

crystallographic  meaning  and  are  defined  in  Appendix  C. 

The  standard  theoretical  and  experimental  procedures 

12  3 

can  be  found  in  detail  elsewhere  '  '  . 


Chapter  I 


THE  HINDERED  ROTOR 


INTRODUCTION 


In  X-ray  diffraction  the  structure  factor  F 

2 

can  be  generally  expressed  as 


hk£ 


pxyz  exPt27T^  h-r]  dx  dy  dz  .  [1] 

x  y  z 

If  the  assumptions  are  made  that  atoms  are  at  rest 
and  the  electron  density  is  zero  elsewhere  in  the  cell, 
the  expression  simplifies  to 

Fhk i  =  ?  fi  exP  [2i\i  (hx±  +  ky±  +  H  z±)  ]  [2] 

where  f .  is  the  scattering  amplitude  of  atom  i  at  an 
angle  0  and  allows  for  the  finite  size  of  the  atom. 

In  general,  atoms  are  not  at  rest,  but  are  in  continual 
motion.  To  allow  for  this,  Debye  derived  a  tempera¬ 
ture  parameter  B  which  was  related  to  the  mean  square 
displacement  of  an  isotropic  harmonic  oscillator.  Then 

Fhk*  became 

Fhk£  =  £  fi  exp[-B  (sine/A  )  2]exp[27ii  (hxi+kyi+£zi)  ]  .  [3] 

In  many  cases  a  better  description  of  the  motion  of  the 

atoms  can  be  obtained  by  assuming  a  general  rectilinear 

4 

harmonic  oscillator  .  This  model  is  conveniently  handled 
in  tensor  form  and  the  six  unique  thermal  parameters  of 
the  general  case  are  related  to  the  root  mean  square 


2 


amplitudes  and  the  direction  cosines  of  the  three 
principal  axes  of  the  ellipsoid  of  vibration.  The 
isotropic  and  anisotropic  models  are  usually  considered 
practical  to  describe  most  crystal  structures. 

The  major  disadvantage  of  these  models  lies  in  the 
large  number  of  parameters  required  to  describe  the 
crystal  structure,  especially  in  large  molecules  such 
as  RhH (CO)  (P  (C^H^) ^ .  Since  the  isotropic  model 
requires  four  parameters  per  atom  (x,y,z,B)  in  the 
general  position  and  the  anisotropic  model  requires 
nine  parameters  per  atom  (x,y,z  and  for  example  3  , 

$22'  ^33'  ^12'  ^13'  ^23)  a  ful1  matrix  refinement  of  a 
large  molecule  would  require  excessive  computer  time 

and  storage.  There  are  two  obvious  solutions  to  this 
problem.  The  first  is  to  use  block  diagonal  refinement^ 
rather  than  full  matrix  refinement.  This  cuts  down 
greatly  on  both  time  and  storage  requirements,  but  has 
the  disadvantage  that  some  large  parameter  correlations 
can  be  ignored,  even  with  an  apparently  sensible  choice 
of  blocks. 

The  second  method  retains  the  use  of  full  matrix 

refinement  and  utilizes  known  geometries  for  certain 

7  8 

rigid  groups  '  ,  the  disadvantage  of  the  method  being 
that  significant  deviations  from  ideal  geometry  will  be 
missed.  In  this  method,  certain  groups  of  atoms  are 


' 


' 

' 

: 

. 
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treated  as  a  unit.  It  is  assumed  that  the  groups 
are  internally  rigid  and  the  temperature  factors  of 
the  atoms  are  dependent  only  on  the  motion  of  the 
rigid  body  as  a  whole.  Such  a  group  may  be  defined 
by  setting  up  an  internal  Cartesian  co-ordinate  system 
within  the  rigid  body.  The  positions  of  individual 
atoms  within  the  group  can  then  be  specified  with 
reference  to  the  internal  co-ordinate  system.  This 
co-ordinate  system  can  be  related  to  the  crystallo¬ 
graphic  co-ordinate  system  by  means  of  a  transformation 
matrix.  Since  the  atoms  within  the  group  are  assumed 
to  have  fixed  positions  relative  to  the  internal  co¬ 
ordinate  system,  it  is  only  necessary  to  refine  the 
origin  and  orientation  of  the  co-ordinate  system  with 
respect  to  the  crystallographic  axes,  as  well  as  refining 
temperature  factors  for  the  individual  atoms  and/or  the 
rigid  body  as  a  whole.  Use  of  such  a  model  greatly 
reduces  the  number  of  parameters  to  be  refined. 

Use  of  isotropic  temperature  parameters  for  the 

atoms  is  based  on  the  assumption  that  the  rigid  body  is 

undergoing  only  relatively  small  oscillations  from  its 

9 

mean  position  .  The  opposite  extreme  occurs  if  the 
rigid  body  is  rotating  freely.  Expressions  for  the 
scattering  of  X-rays  from  such  groups  have  been  calculated 
by  Bijvoet  and  Ketellar^  and  Zachariasen^ . 


4 


Intermediate  between  the  two  extremes  of  free 
rotation  and  no  rotation  is  the  case  where  the  rotation 
of  the  group  is  hindered.  Cruickshank  found  this  to 
be  the  case  in  benzene  at  -3°C.  After  anisotropic 
refinement  of  the  carbon  atoms  ,  he  found  that  the 
tangential  motion  of  the  atoms  was  much  greater  than 
the  radial  motion  indicating  an  oscillation  of  the  benzene 
ring  about  the  hexad  axis.  In  order  to  describe  this 
motion  better,  he  assumed  that  the  benzene  ring  was  a 
rigid  body.  Then  its  motion  could  be  described  in  terms 

of  two  tensors,  T  describing  the  translational  motion, 

13 

and  oo  associated  with  oscillations  .  Using  this  method 
of  analysis,  he  showed  that  the  major  ring  oscillation 
was  about  the  hexad  axis  with  a  RMS  (root  mean  square) 
angular  oscillation  of  7.9°.  From  this  it  was  also 

possible  to  calculate  the  apparent  shortening  of  the  C-C 

14 

bond  due  to  the  motion  of  the  ring  .  The  method  was 
later  extended  by  Schomaker  and  Trueblood  to  the  case 
where  there  was  no  pre-ordained  centre  of  rotation  and 
allows  for  a  coupling  of  rotation  and  translation  via  the 
introduction  of  a  third  tensor  S.  The  Cruickshank  T  and 
oo  tensor  elements  have  been  used  as  model  parameters 

16 

within  a  structure  factor  and  least  squares  programme 
although  strictly  speaking  the  site  symmetry  of  the  group 
is  limited.  The  general  case  of  Schomaker  and  Trueblood 
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has  not  been  incorporated  into  a  structure  factor  and 
least  squares  programme  and  its  application  is  currently 
limited  to  the  analysis  of  the  results  of  an  anisotropic 
independent  atom  refinement. 

A  different  approach  to  the  problem  of  hindered 

.  1 7 

rotation  was  taken  by  Lipscomb  and  King  in  1949,  at 

which  time  the  anisotropic  tensor  description  was  not 

known.  They  derived  an  expression  for  the  calculation 

of  the  intensity  of  X-ray  scattering  from  rigid  groups 

undergoing  hindered  rotation.  The  assumption  was  made 

that  the  rigid  group  rotates  about  a  fixed  axis  with  a 

sinusoidal  hindering  potential.  The  method  was  used 
1 8 

by  Prewitt  in  1965  when  he  treated  the  cyclopenta- 
dienyl  rings  in  [C5H5FeS]4  as  hindered  rotors.  He 
assumed  that  each  cyclopentadienyl  ring  formed  a  regular, 
planar  pentagon  with  D5h  symmetry  and  that  the  motion  of 
the  atoms  is  predominantly  in  the  plane  of  the  ring.  The 
other  motions  of  the  group  were  described  by  a  single 
isotropic  temperature  factor  B,  this  being  a  modification 
of  the  original  King  and  Lipscomb  model.  This  model  was 
refined  by  means  of  two  matrices  in  a  block  diagonal 
refinement;  one  4x4  matrix  for  x,y,z  and  B,  and  a 
2x2  matrix  for  the  radius  of  the  ring,  and  the  relative 
barrier.  Use  of  this  model  gave  cyclopentadienyl  C-C 
bond  lengths  in  good  agreement  with  the  average  C-C  bond 


- 


- 
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length  obtained  when  corrections  for  thermal  motion 
had  been  made  for  the  individual  atom  models.  It 
also  required  far  fewer  parameters  to  describe  the 
system  although  this  resulted  in  a  higher  R  factor; 

6.9%  for  the  hindered  rotor  as  compared  to  5.8%  for 
the  individual  atom.  Also,  the  hindered  rotor  model 
had  the  advantage  that  it  placed  the  electron  density 
on  the  arc  of  a  circle  giving  a  curved  distribution 
rather  than  attempting  to  describe  the  distribution  in 
a  rectilinear  manner  as  in  the  conventional  isotropic 
and  anisotropic  models.  However,  there  were  disad¬ 
vantages.  The  particular  blocking  of  the  least  squares 
matrix  was  unfortunate  in  that  the  correlation  between 
the  two  parameters  associated  with  a  description  of  the 
thermal  motion  was  ignored.  This  caused  severe  over¬ 
shifting  and  necessitated  the  introduction  of  damping 
factors.  Because  the  ring  orientations  were  well 
defined  from  the  previous  individual  atom  refinement, 
these  were  not  refined  and  the  programme  was  developed 
with  a  limited  scope.  However,  the  utility  of  the  model 
had  been  demonstrated  for  a  routine  organometallic 
structure . 

In  order  to  refine  n  parameters  by  means  of  full 
matrix  least  squares  refinement,  a  symmetrical  matrix 
containing  (n/2) (n+1)  unique  elements  is  required.  Both 
the  computing  time  required  and  the  storage  demands  for 


' 
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a  least  squares  cycle  are  proportional2  to  (n/2) (n+1) . 
Thus,  any  decrease  in  the  number  of  parameters  to  be 
refined  will  result  in  a  proportionately  greater  decrease 
in  time  and  storage  requirements  and  in  the  computing 
cost.  It  was  therefore  deemed  more  economical  to 
develop  a  programme  that  allowed  the  ring  orientations 
to  be  refined  rather  than  to  obtain  these  parameters 
from  a  complete  individual  atom  refinement. 

In  the  following  work,  the  rigid  body  expressions 
available  in  the  rigid  body  modification  of  the  basic 
full  matrix  least  squares  programme  SFLS5  by  C.T.  Prewitt 
(1966)  were  used  in  conjunction  with  Lipscomb  and  King's1^ 
hindered  rotor  expressions.  SFLS5  was  then  further 
modified  to  calculate  the  contribution  to  the  structure 
factors  by  a  group  undergoing  hindered  rotation  and  to 
do  a  full  matrix  least  squares  refinement  on  the  input 
parameters . 


THEORETICAL 


complex  scattering  power  of  the  i^1 
atom 

average  complex  scattering  power  of 
atoms  in  a  ring 

probability  function  of  the  ring  at 
an  angle  0 ' 

reciprocal  lattice  vector  corresponding 
to  the  Miller  indices  h,k, £ 
instantaneous  position  vector  of  the 
atom  which  may  be  resolved  into  two 
components : 

vector  from  the  real  lattice  origin 
to  the  centre  of  rotation,  and 
vector  from  the  centre  of  rotation  to 
the  instantaneous  position  of  the 
atom 

potential  barrier  to  rotation 
potential  at  rotation  angle  0 ' 
number  of  minima  in  the  potential 
instantaneous  rotation  angle  measured 
from  the  x  axis  of  the  internal  Cart¬ 
esian  co-ordinate  system  in  the 
rigid  body 


' 

, 

•  ■ 
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D,E 


R 

X(i) 
Y(i)  • 

z(i)  , 

[L] 


q-L  to  q6 


fractional  co-ordinates  of  the 

centre  of  gravity  of  the  i  n  rigid 
body 

rotation  angle  at  which  the  potential 
is  a  minimum 

angles  (in  radians)  defining  the 
orientation  of  the  plane  of  the 
ring.  The  direction  cosines  (with 
respect  to  the  chosen  internal 
Cartesian  co-ordinate  system)  of 
the  plane  are : 

m^  =  cosD  sinE 
m^  =  -sinD 

and  m^  =  cosD  cosE 

o 

radius  of  the  rigid  body  in  A 


o  th 

Cartesian  co-ordinates  in  A  of  i  atom 

transformation  matrix  from  the  internal 
co-ordinate  system  of  the  rigid  body 
to  conventional  crystallographic 
co-ordinates 

non-zero  elements  of  the  matrix  [L] 

* 

q^  =  l./(a  sin3  siny  ) 

* 

=  coty  /b  sina 


10 


B  . 

l 


0 

X 

f  . 
1 

T 

k 

9 

Bd 


* 

cosa  cosy  -cot$ 


q^  =  l./(b  sina) 

=  -c  cota 

q6  =  l./c 
population  factor  of  the  iLIi  ring 
Debye  temperature  factor  for  the  ith 
ring 

conventional  Bragg  angle 

wavelength  of  X-radiation  used 

form  factor  for  atoms  of  the  ith  ring 

absolute  temperature 

Boltzman's  constant 

imaginary  number  =  /^T 

relative  barrier  Bd  =  V  /2kT 


smy 


+-U 


The  rigid  body  expression  of  M.J.  Bennett  and 
B.M.  Foxman  assumes  a  rigid  ring  system  based  on  a 
regular  planar  polygon  of  order  n  possessing  Dnh  symmetry. 
A  right  handed  Cartesian  co-ordinate  system  is  defined 

such  that  the  origin  lies  at  the  centre  of  the  polygon, 

.  * 
with  the  x  and  z  axes  being  parallel  to  the  a  and  c  axes. 

If  the  ring  lies  in  the  XY  plane,  then  the  co-ordinates 


■ 


i  - 

. 
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of  the  first  atom  can  be  defined  as  (R,  0,  0).  The 

th 

position  of  the  i  atom  is  then 

X  (i)  =  R  cos  2 it  (i-l)/n 

Y  (i)  =  R  sin  2tt  (i-l)/n  [4] 

and  Z (i)  =  0 

with  respect  to  the  Cartesian  co-ordinate  system.  The 

direction  cosine  of  the  ring  normal  in  its  standard 

orientation  is  (0,  0,  1),  Rotating  this  polygon  about 

th 

the  z  axis,  the  i  vertex  becomes 

X  (i)  =  R  cos  U  +  2tt  (i-l)/n) 

Y(i)  =  R  sinU  +  2  7T  (i-1)  /n)  [5] 

Z(i)  =  0  . 

The  direction  cosines  of  the  plane  are  still  (0,  0,  1). 

If  this  is  followed  by  a  rotation  of  D  radians  about 

the  x  axis  and  E  radians  about  the  y  axis,  the  direction 

cosines  of  the  plane  become  (cosD  sinE,  -sinD,  cosD  cosE) . 

th 

The  co-ordinates  of  the  i  atom  are  then 

X(i)  =  R[cos(£  +  2u(i-l)/n)  cosE 

+  sin  U  +  2iT(i-l)/n)  sinD  sinE] 

Y  (i)  =  R  [sin  (£  +  2ir(i-l)/n)  cosD]  [6] 

Z  (i)  =  R[-cos(C  +  2ii(i-l)/n)  sinE 

+  sin(£  +  27r(i-l)/n)  cosE  sinD] 
Transformation  to  crystallographic  co-ordinates  using 


12 


matrix  [L]  plus  translational  correction  gives  the 

x, 

fractional  co-ordinates  of  the  1  atom  as 

=  q  ^  R[cos(£  +  2iT(i-l)/n)  cosE 

+  sin  (?  +  2TT(i-l)/n)  sinD  sinE]  +  x 

c 

=  R[cos(£  +  27T(i-l)/n)  cosE 

+  sin(C  +2ff(i-l)/n)  sinD  sinE] 

+  q.  R[sin(£  +2TT(i-l)/n)  cosD]  +  y 
*  c 

=  q^  R[cos(?  +  27r(i-l)/n)  cosE  [7] 

+  sin  (£  +2n(i-l)/n)  sinD  sinE] 

+  q j.  R[sin(?  +2n(i-l)/n)  cosD] 

+  q^  R[-cos(C  +  2u(i-l)/n)  sinE 
+  sin  +  2ir  (i— l)/n)  cosE  sinD]  +  z 

c 

where  q^  are  the  elements  of  the  transformation  matrix 
[L] .  Replacing  the  fixed  angle  £  by  the  continuous 
variable  0  '  ,  the  dot  product,  2Trh*y,  used  in  later 
expressions,  becomes 


2-rrh‘V  =  a  cos9'+  dsin0  ' 

[8] 

where 

o  =  2ttR[  (h  q^  +  k  q^  +  Jl 

q^)cosE  - 

£  q^  sinE] 

[9] 

and  d  =  2ttR[  (h  q^  +  k  q^  + 

q^)sinD  sinE 

+  (k  q^  +  £  q^ )  cosD  +  £  q^  cosE  sinD], 


x  (i) 

y  (i) 

and  z (i) 


!<) 


' 

' 
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If  a'  is  defined  as  a'  =  /u  +  d^ 

and  n  =  arc  cosc/a' 

=  arc  sin  d/a' 

then  the  dot  product  27Th*v  becomes 


[10] 


[11] 


2nh *v  =  a'cos (6 1  -  n ) •  [12] 

These  rigid  body  expressions  can  now  be  used  in  calcu¬ 
lating  the  contribution  of  groups  in  hindered  rotation 

to  the  structure  factor.  The  derivation  parallels 

17 

that  of  Lipscomb  and  King 

The  contribution  of  the  hindered  rotor  groups  to 
the  structure  factor  can  be  expressed  as 


F 

hk£ 


[13] 


l 

the  summation  being  over  all  atoms  in  the  unit  cell  which 
are  members  of  a  group  undergoing  hindered  rotation.  The 
instantaneous  complex  scattering  power  of  an  atom  in  a 
ring  is  given  by 


gi  =  f^  exp  [ 2tt£  h*r]  [14] 

=  f^  exp[27T'£  h  •  k]  exp  [  2iri  h*v]  . 

If  the  probability  distribution  is  given  by  P  (6 ' ) ,  then 

g  =  f^  exp[2nt  h-k]  • 


*  2  7T 

P(e')  exp  [2-ni  h*v]d0'.  [15] 

o 

Using  a  potential  function 

V  =  -h  VQ  cos  n(6’  -  £) 


[16] 


' 

■ 
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and  assuming  a  Maxwell— Boltzman  distribution  for  the 
probability  function,  P  (0 1 )  ,  then 


exp  [V  cos  n(0'  -  £)/2&T] 

P(6-)  =  — - 2 - -  .  [17] 

exp [V  cos  n (9 '  -  i) /2kT]dQ  ' 
o 


Substituting  this  expression  for  P(Q')  into  equation  [15] 
and  using  the  dot  product  expression  of  [12] ,  one  obtains 


g  =  f^  exp[27Tt  h*k]  *  [18] 

■2tt 

exp  [Bdcos  n  (0  '  -  £)  +  ^a'cosCe'  -  n)]d0' 
o _ _ 

"  2  n  _________  _ 

exp  [Bdcos  n(0'  -  E, )  ]  d0  ' 

o 


where  Bd  =  V  /2kT . 

o 

Defining  M^(a',Bd)  = 


■2tt 

exp  [Bdcos  n(0'  -  5)  +  ia  '  cos  (0  '  -  £  )  ]  d0  ' 


•  2  TT 

exp [Bdcos  n(0'  -  5)]d0' 
‘  o 


[19] 


then  g  becomes 


g  = 

(a' , Bd) 
n 

T7  4-  20 

Watson  gives 


f  ^  exp[2iri  h*  k]  M^(a',Bd). 
can  be  expanded  into  an  infinite  series, 
the  following  expression: 


[20] 


i 


' 


. 


. 
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exp[s  cost]  -  I0(s)  +  2  l  I  (s)cos  mi 

m=l 

00 

=  V”(s)cos  mT  t21] 

m=0 

Where  em  =  1  when  m  =  and  e  =  2  when  m  ^  0.  J  (s) 
m  m  m 

is  a  Bessel  function  of  order  m  and  I^(s)  is  a  hyperbolic 

Bessel  function  of  order  m.  Using  these  expressions,  the 

numerator  of  (a' ,Bd)  becomes 

n 


*2  7T 

[Ze:pIp(Bd)cos  p(ne  '  “  ]  )cos  m(0  '  -5)6'] 


=  He  £  I  Bd  I  a' 
m  p  p  m 

pm  ^  * 


2  7T 


[22] 


cos  p (ne '  -  c ) cos  ( 0 '  -  n)de'. 


By  expanding  the  cosine  term  the  numerator  can  be  written 
as 


Z  Zz  z  I  (Bd)  I  (£a *  )  • 
^  m  p  p  '  mv  ' 
pm  c  c 


*2tt 

(cos  pn0 '  cos  pC  cos  me'  cos  mn  [23] 

o 

+  sin  pn0 '  sin  pC  cos  m0 '  cos  mn 

+  cos  pn0 '  cos  pC  sin  mO '  sin  mn 

+  sin  pn0 '  sin  pC  sin  m0 '  sin  mC)d0'  . 

Making  use  of  (1)  summation  of  integral  equalling  the 

integral  of  summation,  and  (2)  the  orthogonality  of  sines 
and  cosines,  this  expression  simplifies  to: 


. 


■ 
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277  Z  epIp(Bd)Ipn(t:a,)  (cos  cos  pnn  +  sin  sin  Pnn )  .  [24] 

P 


Using  the  relationship 

J  (-is)  =  (-i)ml  (s) 
m  m 

2  0 

given  by  Watson  the  numerator  becomes 


[25] 


2n  ^pnI  (Bd)J  (a')cos  (pnU  -  n))  [26] 

p  ir  Jr  xrAA 

and  by  the  same  method  the  denominator  is  2ttI  (Bd)  , 

o 

M^(a',Bd)  then  becomes 

ZepiPnVBd)Jpn(a,)C°S(pnU  -  n)> 

2 - -  •  [27] 

VBd) 

Using  this  expression,  the  contribution  to  the  structure 
factor  from  the  hindered  rotors  can  be  written  as 


Fhlc«  =  Efi9P,  exp[-B  SJ-n  S]  exp  [  2iri  h-k]M  1  (a!,Bd.)  [28] 

ms.  jo.  xx  -L  ,  z  —  —  n.  11 

A  1 

the  summation  being  over  all  the  rings  undergoing  hindered 
rotation . 


Programming 

The  programme  SFLS5  was  modified  to  calculate  the 
contribution  of  a  hindered  rotor  to  the  structure  factor. 
These  modifications  were  mainly  in  the  addition  of  two  new 
subroutines  ROTOR  and  BESSEL,  as  well  as  modifications  to 
the  input  and  output.  It  can  be  seen  that  there  are  two 


' 
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possible  routes  by  which  to  calculate  the  M  function: 

(1)  by  means  of  the  integral  expression,  equation 

[19] ,  or 

(2)  by  use  of  the  infinite  series  summation 

form  (Bessel  functions) ,  equation  [27] . 

The  first  of  these,  the  integral  form,  can  be 

19 

calculated  by  means  of  Weddel's  Rule  .  If  one  divides 
the  range  of  integration  into  7  intervals,  then 


rxo+6h 

j  0(x)dx  =  Ig-[yo  +  5yx  +  Y2  +  6y3  +  Y4  +  5y5  +  yg]  [29] 

X 

o 


This  can  be  expanded  for  use  with  (6n+l)  intervals,  where 

;b 

n  is  an  integer.  An  integral  of  the  form 


f (x)  ,  .  _ 

e  g(x)dx 


can  be  expanded  into  a  series  of  integrals  where  f (x)  and 
g(x)  are  two  functions  of  x. 


b  f(x)  ,  .  ,h 

e  g (x)  dx  = 

a  '  a 


[1  +  f  (x)  +  -jt  (x)  +  ...]  g(x)  dx 


g  (x)  dx  + 


f  (x)  g (x)  dx 


[30] 


+  h 


f  (x)  g (x)  dx  +  . . 


which  hopefully  are  easier  to  integrate.  Thus  the  M^(a',Bd) 
function  can  be  broken  up  into  a  series  of  integrals  which 
can  be  calculated  by  means  of  Weddel's  Rule,  or  integrated 
analytically.  However,  it  was  found  that  even  integrals 


■ 


' 

■ 
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of  the  form  g(x)  f  (x)  dx  became  extremely  compli- 

cated  (Figure  1)  and  required  a  minimum  of  67  intervals. 
Since  this  is  an  extremely  slow  method  of  calculation, 
it  was  used  only  as  a  check  on  the  results  of  the  infinite 
Bessel  function  series  calculations. 

The  infinite  series  form  of  (a ' , Bd)  is 

n 


r  I  (Bd) 

M^(a\Bd)  =  I  e  iP"  J  (a  • )  5^—,  cos  (pn  U  -  n))  [27] 

p  ^  *  o 


where  I  (Bd)  =  Z  y,P- 

p  m  m!  (m  p)  : 


[31] 


and  J  (a ' )  =  Z 
pn 

^  m 


(-1) m (a 1 /2 ) 2m+pn 

m] (m+pn) I 


[32] 


20 

J  (a')  being  a  Bessel  summation  and  I  (Bd)  a  hyperbolic 

P-H  p 

Bessel  summation.  It  can  be  seen  that  whether  any  term  in 

the  summation  is  real  or  imaginary  depends  completely  on 

the  value  of  ipn;  all  other  functions  are  real.  Thus, 

if  n,  the  order  of  the  ring,  is  even,  pn  is  always  even 

and  £pn  is  ±1;  (a' ,Bd)  is  then  a  real  summation.  If  n 

n 

is  odd,  then  pn  alternates  even  and  odd  as  p  increases 
causing  ipn  to  alternate  among  ±1  and  ±i .  M^(a',Bd)  is  a 

complex  function  which  can  be  broken  up  into  real  and  imag¬ 
inary  contributions : 


. 


' 


■ 
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Figure  1 


1  2 

2  g(x)f  (x)  versus  x 


where  g(x)  =  cos [a ' cos  (x-n ) ] 
f (x)  =  Bdcos  n  (x-£ ) 

The  following  values  were  used: 


h  =  1 

k  =  2 

II 

H* 

5  =  0.23 

D  = 

0.60  radians 

n  =  5 

E  = 

0.028  radians 

Bd  =  2.5 

\ 

R  =  1.20A 

=  0.10 

q2  =  0.20 

q^  =  0.10 

q4  =  o.io 

q5  =  0.27 

q6  =  0.10 

■ 

(x)6(xy 
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M^(af,Bd) 


Mr  +  iMj. 


P 

+ 


(  l)pn/2£p  Jpn(a’)^TBdy  cos  pn(S-n)  [33] 


I  (Bd) 


even 


pn-1 

1  (-1)  2  ep  JDn(a')  I  (Bd)  cos  pn(C-n). 

p  odd  v  p  P 


Fhki  comP^ex  anc^  can  also  be  broken  down  into  real  and 

imaginary  terms,  generally  called  A,  .  _  and  B,  ,  . 

hk£  hk  i 


=  Zf. 

•  l 
l 

exp  (  2tt i 

h-k)  (a' 

—  —  n 

<Bd) 

=  Z  f  . 
i  1 

(cos  27Th*k  +  isin 

27Th’k)  (M  +  il\ 
R 

[34] 

• 

•  • 

A,. 

hk£ 

=  zf. 
i 1 

(M  cos 

27Th*k  -  M-j. 

sin 

27Th-k) 

[35] 

and 

Bhk£ 

=  Zf. 

1 

(M  cos 

27ih‘k  +  M 

R 

sin 

27rh  •  k) 

[36] 

Thus 

for  n 

odd  it 

is  necessary  to 

calculate 

both  M 

R 

and 

Mj  as 

well 

as  the  sin  and  cos 

terms  even  for 

the  centro- 

symmetric 

case . 

For  the  special 

case 

of  even  membered 

rings  in  a 

centrosymmetric  space  group 

then 

A,  . 
hk£ 

=  z  f. 
i  1 

M  cos 
K 

2-rrh  •  k 

[37] 

In  a  centric  space  group  with  no  anomalous  dispersion 

Bhk£  =  °*  This  ProPertY  was  used  as  a  further  check  on 
the  calculation  of  and  M  since 

cos  2irh*k  =  -M  sin  2nh-k  . 


[38] 


' 
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In  order  to  do  a  least  squares  refinement,  it  is 
necessary  to  calculate  the  derivatives  of  the  M^(a',Bd) 
function.  The  ring  is  described  by  means  of  ten  para¬ 
meters,  and  analytical  expressions  for  the  derivatives 
of  F^^  with  respect  to  each  parameter  were  obtained. 
Those  of  the  population  factor  gp,  the  co-ordinates  of 

the  centre  of  gravity  x  ,  y  ,  z  and  the  overall  temp- 

c  c  c 

erature  parameter,  B,  were  straightforward,  since  they 
are  not  involved  in  the  calculation  of  the  M  function. 
Since 


Fhk*  =  E  fi  exP 


2„n 


•B.sin  0 
1 


exp[2irf  h*k]ivr  (a'  ,Bd)  [28] 


then  the  derivative  with  respect  to  the  population  factor 
becomes 


3  F 


3  (gp) 


f±  exp 


-B . sin 
i 


2 


exp[27if  h*k]iyr  (a '  ,Bd) 


[39] 


Using  the  expansion  of  h*k  to  (hx  +  ky  +  z  )  and 


c . 

l 


taking  the  derivative  of  F  with  respect  to  x 


3F 


3x 


=  gPi  ^  exp 


-B .  sin  0-1 
1 


c . 
1 


(2t xi  h) 


ci 


exp[27Ti  h*k]M^(a*  ,Bd)  . 


[40] 


Analogous  expressions  can  be  written  for  ^ 


3F 


3F 


3  z 


c . 

l 


and 


. 


. 


. 
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where  2i\i  h  is  replaced  by  2tt i  k  and  2tt i  l  respectively. 
Differentiation  with  respect  to  the  overall  temperature 
factor  gives 


9  F 
9  B 


.  2 
sm  0 


exp 


B .  sin 

l 


2  -r 


f^  exp[27T-£  h*k]M^  (a  '  ,Bd) 


[41] 


The  remaining  five  parameters  are  used  in  the 

calculation  of  M^(a',Bd)  and  the  analytical  expressions 

for  their  derivatives  are  somewhat  more  complicated. 

The  relative  barrier,  Bd,  is  found  only  in  the 

terms  I  (Bd)  and  I  (Bd) .  The  derivative  of  M^(a',Bd) 

^  n 

with  respect  to  Bd  is  given  by 


SMMa',Bd)  1 

n  -  =  — - -  [I„  (Bd)  E  ipnJ  (a') 

I  (Bd)  p  ^ 


9  Bd 


3  1  (Bd) 

cos  pn(5-n)  fBd 


91  (Bd) 

o  „  • pn ,  ,  ,  x 

- -  Z  £  ^  J  (a') 

9  Bd  p  pn 

P 


I  (Bd)  cos  pn  )  ] 
P 


Since 


I  (Bd) 
P 


=  z 

m 


(Bd/2 ) 


2m+p 


ml  (m+p) ! 


then 


91  (Bd) 

,  _ _ 

9  Bd 


=  Z 

m 


2m+p  (Bd/2) 2m+P 
Bd  m.'  (m+p)  l 


[42] 


[31] 


[43] 


I 


■ 


■ 
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That  is,  each  term  in  this  summation  is  times  the 

Bd 

original  Bessel  term.  An  exactly  similar  expression 
occurs  for 


31  (Bd) 
o 

3  Bd 


=  E 


2m  (Bd/2 ) 


2m 


Bd  ,  ,.2 
m  (m .  ) 


[44] 


It  was  for  the  calculation  of  such  modified  Bessel  terms 
that  the  subroutine  BESSEL  was  written  rather  than  using 
the  computer  library  function. 

In  the  expression  for  M^(a',Bd),  only  the  term  J  (a1) 
is  a  function  of  the  radius,  R. 


3 M  (a  '  ,  Bd) 
n 

3  R 


i  ;pn  e 

p 


I  (Bd) 
P 

P  i0(Bd) 


cos  n(?-ri) 


3  J  ( a 1  ) 
pn 

3  R 


[45] 


Since 


J  (a'  ) 
pn 


E 

m 


(-l)m(a'/2) 2m+pn) 


ml (m+pn) ! 


[32] 


then 


3  J  (a') 
pn 

3  R 


^  (-1) m (a 1 /2 ) 2m+pn-1  3(a«)  ,  , 

*  - m^rniin)  ! - hi  (2m+pn) 


E 

m 


( —  1 ) m ( a  "  / 2 ) 2m+P ( 2m+pn ) 

m.'  (m+pn)  l 


[46] 


where  a  "  =  a ' /R. 


) 

. 


r 


this  is  a  modified  Bessel 


8  M 

As  in  the  case  of  — 
summation. 

The  derivative 


? (a' ,Bd) 
n 

3Bd 

8M^(a' ,Bd ) 

dl 


l 

P 


p^ 


pn 


e  J 
P  pn 


(a '  ) 


Ip  (Ed) 
Io(Ed) 


is  simply 

sin  (pn  U-n)  )  • 


[47] 


The  most  complex  expressions  for  the  derivatives 
°f  ari-ses  from  those  with  respect  to  the  orientation 

angles  D  and  E.  In  each  term  of  the  M^(a',Bd)  summation 
dpn ( a 1 )  and  cos(pnU-n))  are  functions  of  both  D  and  E. 
Dealing  first  with  the  derivative  with  respect  to  D 


3MC (a' , Bd) 
n 

3D 


VBd) 

VBd) 


Jpn  (a'  }  sin  pn 


in 

3D 


+  cos  pn  (£-n  ) 


z  (~1) m (a  1 /2 ) 2m+Pn  . 
m! (m+pn) ! 


2m+pn  3a' 
———————  #  ■  —  -  ~ 

a'  3D  . 


[48] 


3a' 

3D 


where 


d_ 

a ' 


3d 

3D 


' 


■ 
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3d 


and  3D  =  27TRt  (h  qjL  +  k  q2  +  £  q3) 


-  (k  q^  +  £  q5)  sinD  +  £ 


cosD  sinE 
q^  cosD  cosE]  . 


The  expressions  for  d,a*  and  r\  are  given  in  equations 
[9],  [10]  and  [11]  respectively. 

The  derivatives  of  the  M^(a',Bd)  function  with 
respect  to  E  are  calculated  in  a  similar  manner  with 
the  result: 


3M  (a*,Bd)  I  (Bd) 

n  —  =  Z  ;Pn  _£ 


3E 


Io(Bd) 


pnJ  (a*)  sin  pn(£-n)  -j^r 
c  pn  ^  9E 


+  cos  pn  U~r] )  £ 


( - 1 ) m ( a 1 /2 ) 2m+Pn 
m!  (m+pn) ! 


2m+pn  3a' 
a'  3E 


[50] 


where 


3a* 
3  E 


-2  R[(h  q1  +  k  q2  + 


£  q^ )  sinE  + 


£  qr  cosE] 
o 


3d 

3E 


2iTR[(h  q^  +  k  q2  +  £  q  )  sinE  cosE 


-  £  q^  sinD  sinE] 


in 

3E 


1 _ 

(d/a 1 )  2 


3a’ 

3E 


[51] 


and 
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The  derivatives  of  the  F  with  respect  to  the  last  five 
parameters  can  be  calculated  from  the  derivatives  of  M 
by  application  of  the  chain  rule  to  equation  [28] . 

Once  written,  this  programme  was  tested  by  several 
different  methods.  The  M  function  itself  was  calculated 
by  two  different  methods: 

(1)  the  infinite  series,  the  test  for  convergence 

being  that  the  final  term  was  less  than  0.01% 
of  the  summation,  and 

(2)  the  integral  method  using  Weddel's  Rule19. 

The  results  were  also  compared  with  those  given  by  Lipscomb 
17 

and  King  and  agreed  to  within  the  four  figures  given. 

The  calculations  of  derivatives  were  also  checked  by 
several  methods: 

(1)  by  calculating  ^  analytically  from  the  infinite 

series  form  of  M,  0  being  one  of  the  variables 
D ,  E ,  £  ,  R ,  or  Bd , 

(2)  by  calculating  M  at  various  values  of  0, 

plotting  the  results  and  obtaining  the 
slope,  and 

(3)  by  means  of  the  numerical  analysis  technique 

2 1 

of  forward  differences 

Once  these  tests  indicated  that  the  programming  was  correct, 
the  hindered  rotor  calculations  were  tested  on  several  rings 
found  in  published  crystal  structures.  These  were  (1) 


' 
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12  22 

benzene  ,  ruthenocene  ,  (3)  dicarbonyl— tt— cyclohepta— 

trienyl-a-pentaf luorophenyl  molybdenum23,  (4)  cyclo- 

pentadienyl  manganese  tricarbonyl  ,  and  (5)  benzene- 

.  2  5 

chromium  tricarbonyl  . 

Diagrams  of  these  molecules  produced  using  the 
published  parameters  of  the  individual  atom  refinements 
are  shown  in  Figures  2  to  6 .  With  the  exception  of 
ruthenocene,  the  drawings  show  each  molecule  projected 
onto  the  plane  of  the  ring.  The  diagrams  were  all 
produced  by  C.K.  Johnson's  ORTEP  programme  and  give 
qualitative  indications  of  librations  about  the  ring  axes. 
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Figure  2 

The  structure  of  ruthenocene 
G.L .  Hardgrove  and  D.H.  Templeton 
Acta  Cryst. ,  12,  28  (1959) 


Figure  3 

The  structure  of  benzene 
E.G.  Cox,  D.W.J.  Cruickshank  and  J.A.S.  Smith 
Proc.  Roy.  Soc. ,  A247 ,  1  (1958) 
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Figure  4 

The  structure  of  cyclopentadienyl  manganese  tricarbonyl 

A.F.  Berndt  and  R.E.  Marsh 
Acta  Cryst. ,  16,  118  (1963) 


Figure  5 

The  structure  of  benzenechromium  tricarbonyl 
M.F.  Bailey  and  L.F.  Dahl 
Inorg.  Chem. ,  4,  1314  (1965) 
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Figure  6 

The  structure  of 

dicarbonyl-7T-cycloheptatrienyl-a-pentaf  luoro- 

phenyl  molybdenum 
M.R.  Churchill  and  T.A.  O'Brien 
J.  Chem.  Soc.  (A) ,  1110,  (1969) 


*•> 

. 
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Results 

The  refined  parameters  for  the  hindered  rotor 

calculations  are  given  in  Table  I.  Parameters  for  other 

atoms  of  the  structure  did  not  change  significantly  from 

the  individual  atom  refinement  values  and  are  not  listed. 

In  each  case  the  results  of  the  individual  atom  refinement 

given  in  the  published  paper  were  used  to  obtain  the  input 

to  the  hindered  rotor  calculation.  The  centre  of  gravity 

and  the  orientation  angles  D,  E  and  B,  were  obtained  from 

the  programme  MMMR  (M.J.  Bennett  and  B.M.  Foxman ,  1967) 

which  uses  the  positions  of  the  atoms  in  the  ring  as  input. 

The  radius  was  calculated  from  the  average  distance  between 

the  centre  of  gravity  and  the  atom  positions.  The  barrier, 

Bd,  was  calculated  from  NMR  results  where  available 
27  2  8 

(benzene  and  ruthenocene  ) ;  otherwise  set  at  Bd  =  2. 

The  isotropic  temperature  parameter  for  the  ring,  B,  was 
set  to  be  slightly  greater  than  that  of  the  metal  to  which 
it  was  attached.  In  all  cases,  the  hydrogen  atoms  were 

o 

placed  1.0A  from  the  carbon  atoms.  The  results  are 
discussed  under  the  headings  of  (1)  refined  radius,  (2) 
relative  barrier,  (3)  temperature  factor  of  the  ring, 

(4)  the  R  factor,  and  (5)  time  and  storage  requirements. 
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Standard  deviations  in  parentheses  refer  to  the  last  digit  listed. 


33 


DISCUSSION 

Radius 

In  the  refinement  of  crystal  structures  by  means  of 
the  thermal  ellipsoid  description,  the  origin  of  an  ellipse 
tends  to  refine  to  the  centre  of  gravity  of  the  local 
electron  density.  For  the  case  of  an  atom  moving  along 
a  radial  track,  this  centre  will  not  be  on  the  track.  If 
the  radial  track  is  the  arc  AQB  of  Figure  7 ,  then  the 
centre  of  gravity  would  lie  at  point  P.  The  radius  of 
the  ring  would  show  an  apparent  shortening  of  the  ring 
by  an  amount  PQ.  The  ratio  of  the  apparent  radius  to 
the  true  radius  would  be 


TT 

p  (e ' )  cose'  de ' 
o 

■  7T 

p (e ' )  de 1 

o 


[52] 


In  this  case  a  sinusoidal  distribution  has  been  used  and 
P  (9  * )  can  be  calculated  at  various  6 '  angles  as  a  function 
of  the  relative  barrier  to  rotation  Bd  and  the  order  of  the 
ring.  The  results  are  plotted  in  Figure  8. 

For  benzene,  with  a  relative  barrier  from  the  hindered 
rotor  refinement  of  1.95,  the  ratio  of  the  apparent  radius 
to  that  of  the  correct  radius  would  be  0.988.  The  radius 


< 


i 


* 
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Figure  7 

Radial  track  of  rotation  group 
where  r  radius  of  ring 

OP  =  apparent  radius 
P(e)  probability  distribution 
along  radial  track 
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Figure  8 


Apparent  fractional  radius  due  to  the 
rotation  of  the  ring  as  a  function  of  the 
barrier  Bd,  for  various  values  of  n. 


Apparent  frac¬ 
tional  radius 


•tt/N 

P  (9  )  cos©  d© 
o  _ _ 

r  tt/N 

P  (0  )  d© 
o 


cos0  P(@)d# 


Bd 
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calculated  by  means  of  the  hindered  rotor  method  was 

o 

1.391A.  Using  this  result  and  the  ratio  of  0.988,  one 

o 

would  expect  an  uncorrected  radius  of  1.374A  in  the 

individual  atom  treatment.  This  is  close  to  the  result 
°  12 

of  1.377A  obtained  by  Cruickshank  .  The  hindered  rotor 

o 

value  of  1.391A  is  in  good  agreement  with  Cruickshank ' s 

°  29  ° 

corrected  value  of  1.392A  and  the  Raman  value  of  1.397A. 

Comparable  agreement  is  found  for  the  case  of  cyclopenta- 

dienyl  manganese  tricarbonyl.  A  relative  barrier  from 

the  hindered  rotor  refinement  of  1.50  gives  a  ratio  of 

o 

0.978.  Using  the  hindered  rotor  C-C  distance  of  1.424A 
and  a  ratio  of  0.978  one  would  predict  an  uncorrected 

o 

individual  atom  refinement  C-C  distance  of  1.394A.  Berndt 
and  Marsh24  in  fact  obtained  an  average  C-C  distance  of 
1.394A. 

A  similar  pattern  of  C-C  bond  lengthening  for  the  other 
structures  can  be  seen  in  Table  II.  In  each  case,  the 
C-C  distance  obtained  from  the  hindered  rotor  refinement 
was  multiplied  by  the  appropriate  ratio  to  obtain  a  pre¬ 
dicted  uncorrected  individual  C-C  distance.  As  can  be 
seen  from  Table  II,  these  predicted  values  agree  satis¬ 
factorily  wi  til  the  results  actually  obtained  from  the  indiv¬ 
idual  atom  refinement.  It  would  appear  that  the  hindered 
rotor  model  correctly  allows  for  the  apparent  bond  short- 
due  to  the  motion  of  the  ring,  and  makes  the  riding 


ening 


. 

_ 
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Comparison  of  C-C  bonds 
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correction  unnecessary,  assuming  only  one  dominant 
libration. 


Barrier  to  Rotation 

The  potential  function  for  the  hindered  rotor  was 


V  =  -  i  V  cos  n(e'-£)  [16] 

2  o 

The  barrier  to  rotation  in  kilocalories  per  mole  is  V  . 
The  refined  relative  barrier,  Bd ,  is  related  to  the 
barrier  of  rotation  VQ  by 


Bd 


V  x  1000 
_ o _ 

2  x  1.987  x  T 


[19] 


The  value  of  Bd  is  an  indication  of  the  amount  of 
rotation  the  group  undergoes.  A  barrier  of  infinity 
would  give  no  rotation  while  a  barrier  of  zero  would 
result  in  free  rotation  of  the  rigid  group.  Intermediate 
values  are  directly  related  to  the  root  mean  sguare 
angular  displacement  of  the  group. 


RMS 


f"/N 

0'* 2  P  (0  '  )  d0  ' 
O  _ 

■tt/N 

P (0 ' )  d0 ' 
o 


[53] 


Since  P  (0  '  )  is  a  known  function  of  Bd  (see  equation  [17]), 
angular  displacement  versus  relative  barrier  can  be  made 
for  various  values  of  n  —  the  order  of  the  ring.  Figure  9 


. 


■  r  *  J 
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Figure  9 


Root  mean  square  angular  displacement  versus  relative 
barrier  for  various  values  of  n 
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is  such  a  plot  for  five  up  to  eight  membered  rings. 

Thus  for  Mn(C0)3Cp,  which  has  a  barrier  (Bd)  of  1.49, 

one  would  expect  an  RMS  angular  displacement  of  12° . 

24 

This  is  also  the  value  obtained  by  Berndt  and  Marsh 

13 

by  means  of  Cruickshank ' s  w-T  analysis.  In  the  case 
of  benzene,  the  RMS  angular  displacement  obtained  from 
the  graph  is  8.5°.  This  agrees  satisfactorily  with 
Cruickshank' s  value  of  7.9°12.  RMS  angular  displace¬ 
ments  for  the  other  molecules  are  given  in  Table  III. 

The  values  of  the  relative  barrier  obtained  can  also  be 

27 

compared  to  the  NMR  results  for  benzene  and  rutheno- 
cene28.  Andrews  and  Eades  obtained  an  activation  energy 
for  reorientation  of  the  benzene  ring  about  the  hexad  axis 
of  3.7  ±  0.2  kcal/mole  below  ~240°k.  Unfortunately,  above 
this  temperature  the  barrier  height  drops  rapidly  and  one 
can  at  best  make  a  rough  approximation  to  the  barrier 
from  the  slope  assuming  rotation  is  still  the  dominant 
motion.  At  -3°C,  the  temperature  at  which  the  crystal 
structure  was  done,  this  was  estimated  to  be  1.8  kcal/mole 
giving  a  Bd  of  1.7.  This  is  in  the  same  range  of  values 
as  that  of  the  hindered  rotor  value  of  1.9(1).  The  agree 
ment  of  the  NMR  value  with  that  of  the  hindered  rotor 
appears  to  be  fortuitous.  For  TrCgH^Cr  (CO)  ^  , 

7rC5H5Mn(CO)  3  and  ttC^Mo  (CO)  2CgH5 ,  NMR  results  were  not 
available  with  which  to  compare  the  hindered  rotor  results 


« 


. 


' 

■ 
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The  barriers  (Bd)  obtained  for  these  compounds  were  of 

the  same  order  of  magnitude  as  that  of  benzene  which  would 

be  expected  for  oscillations  of  the  same  order  as  benzene. 

The  results  obtained  are  shown  in  Table  III.  The  value 

obtained  for  ruthenocene  was  greater  than  those  for  the 

other  compounds  and  had  a  much  larger  error  associated 

with  it:  5.0  ±  2.6.  This  did  not  agree  very  well  with 

2  8 

the  NMR  result  obtained  by  Holm  and  Ibers  of 

V  =  2.3  ±  0.2  kcal/moles  which  is  equivalent  to  Bd  =  2.01. 
o 

The  reason  for  this  discrepancy  may  be  in  the  nature  of  the 
potential  function  used.  In  a  cosine  function  the  acti¬ 
vation  energy  required  for  rotation  of  the  ring  should 
equal  V  .  If,  for  example,  the  potential  function  were 
of  the  form  shown  in  Figure  10  which  has  a  cosine  function 
in  the  immediate  area  of  the  potential  minimum  and  a  fixed 
potential  elsewhere,  the  activation  energy  would  be  VQ* 
while  the  "Bd"  described  by  a  purely  cosine  function  would 

be  Bd  =  V  /2kT  and  not  V  '/2/cT.  A  plot  of  P  (0  '  )  versus 
o  o 

6 '  (Figure  11)  for  various  values  of  Bd  shows  that  for 

Bd  =  4.0  or  greater  the  atoms  are  for  the  most  part  confined 

to  within  a  few  degrees  of  the  minimum. 

In  this  case,  Bd  would  be  insensitive  to  the* shape  and 
height  of  the  potential  function  outside  this  region.  Bd 
is  then  a  parameter  describing  the  shape  of  the  function 
close  to  the  minimum  rather  than  a  genuine  measure  of  an 
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A  possible 


P  (0  1  )  versus  6 


Figure  10 


form  of  the  potential  function 


Figure  11 


'  for  a  cosine  distribution  function 
Bd  =  4.0 
Bd  =  2.0 


-  Bd  =  1.0 


Bd  =  0.5 


080 


a. 


8  (DEG.) 
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activation  energy. 

Temperature  Factor  B 

In  the  case  of  hindered  rotors,  the  assumption  is 

being  made  that  most  of  the  motion  of  the  ring  is  in  the 

plane  of  the  ring  and  tangential  to  it.  In  the  case 

where  the  ring  is  strongly  bonded  to  a  metal  atom,  one 

might  expect  the  ring  to  ride  on  the  metal.  If  this  were 

the  case,  the  temperature  factors  for  the  rings  as  a  whole, 

B,  should  be  very  similar  to  the  isotropic  B  of  the  metal 

atom,  B_. .  This  indeed  seems  to  be  the  case  as  can  be  seen 
M 

in  Table  IV.  In  each  case,  the  temperature  factor  seems 
to  be  slightly  greater  than  that  of  the  metal  to  which  it 
is  attached.  This  difference  could  result  from  out  of 
plane  oscillations  by  the  ring  which  are  not  allowed  for 
in  the  hindered  rotor  motion,  and  also  by  librations  of 
the  molecule  as  a  unit. 

The  observation  that  the  isotropic  temperature  factor 
of  the  ring  is  slightly  greater  than  that  of  the  metal  can 
be  utilized  to  obtain  a  starting  value  of  the  barrier. 

This  is  desirable  as  it  is  not  always  possible  to  calcu¬ 
late  a  relative  barrier,  Bd,  from  the  NMR  result,  nor 
would  the  results  be  necessarily  meaningful. 

Application  of  the  hindered  rotor  model  would  gener¬ 
ally  be  preceded  by  individual  atom  isotropic  refinement 


" 
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or  by  rigid  body  refinement.  Both  of  these  methods 
would  give  and  an  average  temperature  factor  for 
atoms  in  the  ring.  B.  .  These  results  can  be  utilized 
to  obtain  starting  values  for  the  hindered  rotor  refine¬ 
ment  as  follows: 

If  one  takes  the  case  where  the  major  motion  of  an 
atom  is  tangential  to  the  ring  and  this  direction  is 
parallel  to  the  z  axis,  the  x  and  z  axes  in  an  orthogonal 
co-ordinate  system  being  in  the  plane  of  the  ring,  then 
the  elements  of  the  less  common  tensor  description  for 
thermal  motion  become 
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The  equivalent  isotropic  B^sq  of  the  atom  is 


[54] 


_  _  .  .  ,  .  8it  -2 

B .  =  B  +  0.6  +  — 5 —  u 

iso  M  3  a) 


[55] 


The  root  mean  square  displacement  (in  A)  is 


TT  /N 


6' 2P  (0  '  )  d0  ' 


u  =  R 

03 
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IT  /N 


P  (0  '  )  d0  ' 


o 


[56] 


Substituting  for  uw  in  equation  [55] ,  one  obtains 
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Table  IV 

Comparison  of  temperature  factors  for  ring 

and  metal  atom 


O  2 

B  (A) 

°  2 

bm  (a> 

°  2 

AB  (A) 

7T-C5H5Mn(CO)  3 

3.64  (17) 

3.18 

0.46 

TT-C.H^Cr  (CO)  _ 

DO  3 

3.09  (9) 

2.63 

0.46 

ruthenocene (1) 

2.11(30) 

1.50 

0.61 

ruthenocene (2) 

2.48(30) 

0.98 

tt-C7H7Mo  (CO)  2- 

3.13(13) 

2.60 

0.53 

"C6F5 

Average 

0.60 

[57] 
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Since  P  (0  *  ) 


tt/N 


8tt 


,2 


o 


tt/N 


is  known  as  a  function  of  Bd, 

P  (0  '  )  de  * 

-  can  be  plotted  as  a  function 


p  (e  ' )  de ' 


'  o 


of  the  barrier  Bd  (Figure  12)  for  various  values  of  n. 

o 

In  Figure  13  a  C-C  bond  distance  of  1.392A  is  assumed 
and 


AB  =  B. 

ISO 


(Bm  +  0.6) 


■tt/N 

p  (0 1 )  e '  de 1 

'  o 
tt/N 

P  (0  '  )  de  1 
o 
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is  plotted  against  the  barrier  Bd.  Thus,  knowing  the 

average  B^sq  of  the  atoms  in  the  ring  and  the  isotropic 

temperature  factor  of  the  metal,  AB,  can  be  calculated 

and  a  starting  value  of  the  barrier  obtained  from  Figure 

12  or  13.  For  example,  in  tt -C^H^Cr  (CO)  the  individual 

atom  refinement  gives  B. 

^  iso 


of  4.5  for  carbon  atoms  in  the 
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2 

AB/R  as  a 


where  AB 


Figure  12 


function  of  barrier,  Bd 
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BARRIER  (Bd) 
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Figure  13 


AB  as  a  function  of  barrier,  Bd 

o 

(C-C  distance  of  1.392A  is  assumed) 


2.6 


benzene  ring  and  B. 

iso 


AB 


B.  “  (B*, 

iso  M 


4.5  -  (2.6 


+ 

+ 


0.6) 

0.6) 


1.3 


From  Figure  13  one  would  obtain  a  starting  value 

of  Bd  =  1.7.  The  final  refined  value  of  the  barrier 

was  1.97.  For  iT-C^H^Mn  (CO)  ^  a  starting  value  of  1.5 

would  be  calculated  while  that  for  tt-C_H_Mo  (CO)  _C,_F.- 

/  7  2  6  5 

would  be  1.8.  The  final  refined  values  were  1.50  and 
1.94  respectively.  Although  not  too  accurate  a  method, 
this  has  the  practical  application  that  one  is  able  to 
obtain  a  starting  value  for  the  relative  barrier  that  is 
acceptably  close  to  the  final  result.  This  reduces  the 
number  of  cycles  required  for  refinement,  since  the 
method  of  non-linear  least  squares  requires  a  good  starting 
value.  A  starting  value  of  the  temperature  factor  of  the 
ring  would  be  B  =  B^  +  0.6.  The  use  of  reasonably  correct 
starting  values  of  the  temperature  factor  is  important 
since  there  is  a  high  correlation  between  the  two  para¬ 
meters,  Bd  and  B,  typically  being  of  the  order  of  0.60. 

Thus  a  poor  starting  value  of  one  could  greatly  affect 
the  other  and  slow  down  refinement.  Such  a  high  cor¬ 
relation  between  the  two  parameters  could  explain  the 

1 8 

oscillation  problems  experienced  by  Prewitt  etal. 


as  they  refined  using  block  diagonal  least  squares,  the 
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barrier  and  ring  temperature  factor  being  in  separate 
matrices . 

Test  on  -CF^,  Groups 

The  hindered  rotor  model  was  also  tested  out  on  the 

fluorines  of  the  four  -CF  groups  in  Cs [ Y (HFA)  ] 26 , 

j  4 

where  HFA  represents  the  hexaf luoroacetylacetonate  ion, 

CF3COCHOCF3.  -CF 3  groups  as  well  as  such  groups  as 

perchlorates  are  noted  for  undergoing  large  thermal 

motions.  The  assumption  was  made  that  the  fluorine 

groups  had  D  symmetry  and  that  the  dominant  motion  was 

3h 

in  the  plane  of  the  fluorine  atoms.  The  results  of  the 
hindered  rotor  calculations  are  given  in  Table  V.  Although 
the  R  factor  is  not  as  good  as  that  of  the  individual  atom 
model,  the  hindered  rotor  model  does  lengthen  the  C-F  and 
C-CF 3  bonds  and  brings  them  close  to  the  corrected  values 
obtained  by  Cruickshank,  Jones  and  Walker  in  ammonium  tri- 
f luoroacetate^ .  Since  -CF3  groups  are  less  rigid,  the 
assumption  of  D3h  symmetry  of  the  fluorine  groups  is 
probably  not  very  good  and  is  the  probable  cause  of  the 
poor  structure  factor  amplitude  agreement  of  the  hindered 
rotor  model.  This  view  gains  some  support  from  the 
deformations  observed  in  the  individual  atom  refinement. 


. 
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Reliability  Factor 

Another  test  of  the  validity  of  the  hindered  rotor 
model  is  on  purely  statistical  grounds.  An  increased 
number  of  parameters  being  refined  will  normally  result 
in  a  lower  R  factor.  Going  from  the  hindered  rotor 
model  to  the  individual  atom  model  requires  more  para¬ 
meters  and  should  result  in  a  drop  in  the  R  factor. 

Whether  such  a  drop  is  significant  can  be  decided  by 

comparing  the  ratio  of  R  factors  with  those  given  by 
33 

Hamilton  .  If  the  ratio  is  greater  than  that  given 

by  Hamilton  then  the  increase  in  the  number  of  parameters 

is  statistically  justified  at  the  significance  level 

chosen  (in  this  case  90%).  Thus  the  hindered  rotor  model, 

which  has  fewer  parameters,  would  be  statistically 

RHR 

justified  only  if  - -  <  Hamilton  ratio.  This  appears 

IA 

to  be  true  for  only  u-C^H^Mn (CO)  and  Ru(C^H^)2. 

Time  and  Storage  Requirements 

One  practical  advantage  of  the  hindered  rotor  model 

over  that  of  the  isotropic  or  anisotropic  atom  model  is 

that  it  requires  fewer  parameters.  Since  the  matrix 

size  =  n(n+l)/2  where  n  is  the  number  of  parameters 

refined,  the  storage  requirements  increase  approximately 

2 

in  proportion  to  n  .  Furthermore,  the  time  required  to 
build  the  matrix  and  later  to  invert  it  is  also  related 
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to  the  matrix  size  and  also  increases  as  n^.  This 
effect  is  especially  noticeable  in  large  structures 
such  as  (C5H5CoCO)2(GeCl2)2Fe(CO)  434.  In  this 
structure  refinement  of  both  molecules  in  the  asymmetric 
unit,  with  the  heavy  atoms  anisotropic  and  all  others 
isotropic,  requires  339  parameters,  350K  bytes  of 
storage  and  two  hours  of  computing  time  per  cycle. 

Use  of  the  hindered  rotor  for  the  four  cyclopentadienyl 
rings  drops  the  number  of  parameters  to  295  and  the 
storage  to  300K  bytes.  The  time  per  cycle  is  cut  in 
half  to  63  minutes.  Thus,  although  the  addition  of 
subroutines  to  handle  the  hindered  rotor  calculations 
lengthens  the  programme  (~12K  bytes),  it  is  far  out¬ 
weighed  by  the  core  saved  due  to  the  large  drop  in  the 
number  of  parameters  to  be  refined,  and  concomitant 
reduction  in  storage  assigned  to  the  matrix. 

For  crystal  structures  such  as  benzene  in  which  only 

a  few  parameters  are  refined,  a  large  percentage  of  the 

computing  time  is  spent  in  calculating  structure  factors. 

Because  calculation  of  the  M  function  is  slower  than  the 

corresponding  individual  atom  expressions,  the  cycle  time 

will  be  slower  for  the  hindered  rotor  model  than  for  the 

individual  atom  model.  However,  for  larger  structures 

such  as  tt-C7H7Mo(CO)  2C6F5,  in  the  order  of  75%  of  computing 
2 

time  is  spent  in  building  and  inverting  the  matrix.  In 
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this  case,  the  slow  calculation  of  the  M  function  is 

outweighed  by  the  reduction  in  matrix  building  time  and 

the  hindered  rotor  model  is  faster.  The  break  even  point 

appears  to  be  for  structures  of  the  size  of  7r-C._H._Mn  (CO)  - 

b  b  3 

for  which  approximately  the  same  time  was  required  in  both 
me thods . 

Conclusions 

The  development  of  the  hindered  rotor  model  within  the 

framework  of  a  conventional  structure  factor  and  least 

squares  programme  was  undertaken  for  specific  application 

to  7T-C  H  groups  bonded  to  metal  atoms.  The  model  is 
n  n  r 

attractive  in  that  it  emphasizes  the  physically  real 
phenomenon  of  hindered  rotation  of  these  groups  in  the 
solid  state.  The  fixed  form  of  the  potential  function 
does  not  appear  to  be  a  serious . problem  as  long  as  no 
physical  significance  is  attached  to  the  height  of  the 
barrier  to  rotation  that  is  obtained.  Bond  lengths 
derived  from  the  hindered  rotor  model  agree  well  with 
those  obtained  from  independent  atom  models  after  cor¬ 
rections  for  thermal  motion  have  been  made.  The  model  has 
been  found  to  be  competitive  with  independent  atom  models 
as  judged  by  a  statistical  comparison  of  conventional 
crystallographic  R  factors. 
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Chapter  II 


THE  CRYSTAL  STRUCTURE  OF  FERROCINIUM  PICRATE 
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INTRODUCTION 

Once  the  hindered  rotor  programme  was  written 

and  tested  it  was  decided  to  test  the  programme  for 

use  in  a  routine  crystal  structure.  The  compounds 

ferrocene  and  ferrocinium  picrate  both  contain  cyclo- 

pentadienyl  rings.  The  crystal  structure  of  ferrocene 

by  Dunitz,  Orgel  and  Rich^  shows  a  libration  of  the 

cyclopentadienyl  rings  about  the  five  fold  axis,  making 

it  ideal  for  the  hindered  rotor  treatment.  Such 

libration  could  also  occur  in  the  ferrocinium  ion  of 

such  salts  as  ferrocinium  picrate.  Since  use  of  the 

hindered  rotor  model  gives  good  average  bond  lengths, 

the  structural  results  could  be  used  as  evidence  in 

determining  the  bonding  in  such  structures.  For  these 

reasons,  the  data  for  ferrocene  were  re-collected  by 

36 

Kathleen  Simpson  and  the  structure  of  ferrocinium 


picrate  undertaken. 
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EXPERIMENTAL 

The  compound  ferrocinium  picrate  [Fe (C^H.- ) _+] 

5  5  2 

—  o  7 

[P^C^N^O^  ]  was  prepared  by  Ester  Lorah  from 
ferrocene  dissolved  in  a  benzene/ethanol  solution. 
This  was  oxidized  with  quinone  in  the  presence  of  a 
stoichiometric  amount  of  picric  acid.  Her  analysis 
showed 


calculated 


C  H  N 

46.36%  2.92%  10.15% 


found 


46.00%  3.06%  10.40% 


A  thermogram  run  to  150°  -  160°  gave  no  weight  loss, 
indicating  no  solvent  and  no  excess  ferrocene  present 
in  the  crystals.  The  crystals  were  long  black  needles 
extinguishing  along  the  needle  axis.  The  symmetry  of 
rotation,  Weissenberg  (hkO  -  hk3)  and  precession  photo¬ 
graphs  (0k£)  showed  the  compound  to  be  orthorhombic. 
Systematic  absences  of  hkJi :  h+k  =  2n+l  and  hOA:  H  =  2n+l 
allowed  three  possible  space  groups  -  Cmc2-^,  Cmcm  or 
C2cm,  a  non-standard  setting  of  Ama2.  Of  these,  only 

o 

Cmcm  is  centrosymmetric .  The  unit  cell  a  =  12.513 (7)A, 
b  =  20.267 (9 ) A  and  c  =  6.903(7)A  (V  =  1750. 6A3)  was 
obtained  from  least  squares  refinement  of  20  values 
of  carefully  centred,  high  angle  reflections  measured 
on  the  Picker  four  circle  manual  diffractometer.  The 
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density  was  determined  by  flotation,  using  a  mixture 
of  benzene  and  dibromomethane ,  to  be  1.60  gm/cc.  This 
is  consistent  with  a  calculated  density  of  1.57  gm/cc 
for  4  molecules  per  unit  cell.  Since  both  space 
groups  Cmc2^  and  Ama2  have  eight  general  positions 
and  Cmcm  has  sixteen,  this  would  require  the  ionic 
species  to  sit  on  special  positions  such  as  a  mirror 
or  two  fold  axis.  Each  space  group  has  one  or  more 
special  positions.  Cmcm  has  site  symmetries  of  mm 
and  2/m.  Site  symmetry  of  mm  is  possible  for  both 
the  ferrocinium  and  picrate  ions.  It  would  require 
the  ferrocinium  ion  to  be  in  the  eclipsed  formation 
while  the  picrate  would  be  required  to  be  planar. 

Site  symmetry  2/m  is  possible  for  the  ferrocinium 
ion  requiring  a  staggered  arrangement  of  the  cyclo- 
pentadienyl  rings,  but  is  impossible  for  the  picrate 
ion  unless  highly  disordered.  The  space  groups 
Cmc2^  and  C2cm  both  have  site  symmetry  m.  Both 
the  ferrocinium  and  picrate  ions  can  possess  this 
symmetry  as  well  as  the  two-fold  symmetry  found  in 
C2cm.  The  possibilities  of  site  symmetries  in  the 
three  space  groups,  assuming  no  major  disorder,  are 
summarized  in  Table  VII. 

Data  were  collected  using  a  crystal  (0.23  mm  x 
0.057  mm  x  0.062  mm)  mounted  along  the  needle  axis  (c) 
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Table  VII 

Possible  site  symmetry 


Space  Group 

Site  Symmetry 
of  FeCp^+ 

Site  Symmetry 
of  Picrate  Ion 

1 

Cmcm 

mm 

mm 

2 

Cm  cm 

2/m 

mm 

3 

Cmc2^ 

m 

m 

4 

C2cm 

m 

m 

5 

C2cm 

m 

2 

6 

C2cm 

2 

m 

7 

C2cm 

2 

2 

• 
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and  with  the  c*  axis  co-incident  with  the  diffracto¬ 
meter  0  axis.  A  moving  crystal,  moving  counter 
technique  (9/20)  was  used  with  a  scan  width  of  2° 
at  two  degrees  per  minute.  CuK  X-radiation  was 


a 


used,  mono chroma ted  by  an  oriented  graphite  crystal 
(002  reflection) .  The  diffracted  beam  aperture 
used  was  1.0  mm  and  the  pulse  height  analyser  was 
centred  on  the  CuKa  peak  to  accept  90%  of  the 
radiation.  Stationary  background  counts  of  30 
seconds  were  taken  on  each  side  of  the  scan.  Six 
hundred  and  forty-two  unique  reflections  (hk£;  h,  k 
and  Z  positive  only)  were  collected,  of  which  276 
were  rejected  using  the  criteria  (1)  1^0  and  (2) 

1^ 3a  where  I  =  [P-t(B1  +  B2)]  and  a  =  [P  +  t2 (B1  +  B2)]^ 

where  P  peak  count 


background  counts 


and  t  ratio  of  total  peak  count  to  total 


background  count 


Seven  well  distributed  standard  reflections  measured 
periodically  throughout  the  data  collection  showed  no 
decomposition  of  the  crystal.  No  absorption  correction 


was  made  (y(Cu)  =  75.1  cm  ^) . 
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SOLUTION  AND  REFINEMENT 

A  three  dimensional  Patterson  map  gave  the  iron 
position  as  x  =  0.0  and  y  =  0.185,  the  z  co-ordinate 
being  undefined.  Although  simple  interpretation  of 
the  Patterson  does  not  determine  the  space  group, 
site  symmetry  considerations  given  in  Table  VII 
reduced  the  possible  arrangements  from  seven  to  two. 

The  Patterson  solution  eliminates  possibilities  2,  4, 

5,  6  and  7,  leaving  only  1  and  3  for  serious  consid¬ 
eration  and  reduces  the  space  group  choice  to  either 
Cmcm  or  Cmc2^.  The  space  group  Cmc2-^  was  chosen  for 
initial  refinement  because  the  crystal  showed  mm2 
symmetry.  Since  this  space  group  is  polar,  the  z 
co-ordinate  was  fixed  at  z  =  0 . 0 .  A  difference  map 
(Fq  “  Fc,  Fc  phased  on  iron  only)  revealed  the  location 
of  all  other  non-hvdrogen  atoms.  With  the  space 
group  chosen  as  Cmc2^,  the  iron  only  phased  Fourier  map 
is  centrosymmetric  and  has  the  symmetry  of  Cmcm.  If 
the  non-centrosymmetric  space  group  is  correct,  the 
effect  of  the  extra  symmetry  is  to  show  the  true  electron 
density  distribution,  together  with  its  superimposed 
mirror  image.  The  picrate  atoms,  with  the  exception 
of  the  ortho  nitro  group,  were  found  on  the  mirror  plane 
(z  =  0.0).  The  appearance  of  diffuse  ortho  nitro  group 
peaks  on  either  side  of  the  plane  could  result  from 
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either  non-centric  space  group  with  its  superimposed 
mirror  image,  or  centrosymmetric  space  group  with 
disorder  of  the  nitro  group.  The  cyclopentadienyl 
rings  of  the  ferrocinium  cation  were  revealed  as 
smears  of  electron  density  and  possessing  mm  symmetry. 
This  could  also  be  due  to  superposition  of  the  mirror 
image  in  the  non-centrosymmetry  space  group,  or  if 
centrosymmetric,  due  to  a  large  oscillation  of  the 
rings.  Calculation  of  the  structure  factors  using 
all  non-hydrogen  atoms  gave  an  R  factor  of  R1  =  26.7%, 
R2  =  28.3%,  where 


Z  |F 

*  i 


and 


The  weights  used  were  obtained  from  peak  and  background 

2  7 

counts  (weight  w  =  1/a  (F) ,  a(F)  =  (P  +  t  (B^  +  B2) 

+  p2I)\  p  =  0.03)  . 

Four  cycles  of  refinement,  minimizing  the  function 
2 

Zw ( | F  | — | F  |)  ,  with  all  atoms  isotropic  and  the  cyclo- 

o  c 

pentadienyl  ring  as  a  hindered  rotor  brought  the  R 
factor  down  to  R^  =  11.6%. 

The  final  orientation  of  the  cyclopentadienyl  ring 
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was  18°  from  the  orientation  required  by  the  additional 
mirror  of  Cmcm.  There  remained  the  distinct  possibility 
of  refinement  into  a  false  minimum  corresponding  to  the 
pseudo-mirror  at  z  =  0.0.  Also,  in  a  non-polar  space 

o  o 

group  such  as  Cmc2  ,  anomalous  dispersion  (in  this 
case  for  the  iron  atom)  can  produce  significant  co¬ 
ordinate  changes.  Four  models  were  then  refined  with 
the  inclusion  of  anomalous  scattering  for  iron: 

(1)  Previous  model;  refinement  terminated  with 

Bd  negative. 

(2)  Mirror  image  of  previous  model;  R  =  11.9%. 

(3)  Ferrocinium  ion  as  in  (1) ,  picrate  as  in  (2) ; 

R  =  11.8%. 

(4)  Mirror  image  of  (3) ;  refinement  terminated 

with  Bd  negative. 

Refinement  of  models  (1)  and  (4)  resulted  in  negative 
values  for  Bd  indicating  that  the  mirror  image  of  the 
cyclopen tadienyl  ring  was  preferred.  This  suggests  that 
the  omission  of  anomalous  scattering  favoured  the  choice 
of  a  false  minimum  for  the  cyclopentadienyl  ring  in  the 
initial  model. 

For  the  ferrocinium  ion,  it  is  reasonable  to  expect 
the  centres  of  gravity  of  the  rings  and  the  iron  atom  to 
be  approximately  collinear.  In  this  crystal  structure, 
this  would  require  the  z  co-ordinate  of  the  ring  centre 
to  be  zero.  However,  models  (2)  and  (3)  showed  a  diff- 
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o 

erence  in  z  corresponding  to  -0.03A. 

•  3  R 

Cruickshank  and  MacDonald's  formula 


A z  (in  A) 


2 


ttS 

max 


max 


where 


S 


2sin0 

A 


predicts  a  difference  in  the  z  co-ordinate  of  0.03A  in 
the  incorrect  model  which  agrees  with  the  observed  shifts. 
This  suggests  that  models  (2)  and  (3)  are  also  wrong  and 
that  none  of  the  cases  considered  so  far  are  correct. 

The  structure  was  then  refined  using  the  centro- 
symmetric  space  group  Cmcm.  This  group  requires  a 
mirror  plane  through  the  plane  of  the  picrate  ion 
and  the  iron  atom.  Isotropic  refinement  of  all  atoms 
brought  the  R  factor  to  =  10.7%,  R2  =  13.7%,  with 
disorder  allowed  for  the  ortho  nitro  group.  These 
atoms  (Nl,  02  and  03)  were  given  half  weight  above  and 
below  the  mirror  plane  for  the  following  reasons: 

(1)  The  initial  difference  map  showed  considerable 
smearing  of  the  peaks  for  these  atoms  off  the  mirror 
plane. 

(2)  When  confined  to  the  mirror  plane,  the  B_ _ 
values  obtained  were  very  large,  suggesting  that  the 
atoms  were  in  fact  off  the  z  =  \  plane,  and  that  the 
ortho  nitro  group  was  not  planar  with  the  rest  of  the 
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picrate  group. 

(3)  If  01,  02  and  03  were  all  in  the  plane, 
there  would  be  a  very  short  01-02  contact  of  2.64A, 
and  an  even  shorter  02-03  contact  of  1.90A,  which  is 
much  shorter  than  the  Van  der  Waal's  radius  sum.  These 
short  contacts  would  be  lengthened  if  02  and  03  were 
off  the  mirror  plane. 

3  9 

(4)  Work  by  Maartmann-Moe  on  the  structure 
of  potassium  and  ammonium  picrates  showed  that  the 
ortho  nitro  groups  are  twisted  25°  -  27°  out  of  the 
plane  of  the  benzene  ring,  making  a  non-planar  picrate 
ion  a  plausible  structure. 

(5)  Study  of  the  packing  diagram  (Figure  15) 
suggests  that  the  nitro  group  could  be  accommodated  in 
either  the  "up"  (z(03)  >0.25)  or  "down"  position. 

In  the  "up"  position,  03  is  within  Van  der  Waal's 
contact  of  H14  of  the  neighbouring  ferrocinium  ion, 
while  in  the  "down"  it  is  the  same  distance  from  H13 . 
Similarly,  02  is  equidistant  from  Hll  and  also  01  in 
either  position.  This  would  suggest  that  the  picrate 
ion  crystallizes  with  a  random  arrangement  of  "up" 
or  "down"  meta  nitro  groups.  The  disordered  centric 
model  was  chosen  rather  than  a  non-centric  model  as: 

(1)  no  clear  distinction  could  be  made  among 
the  four  possible  arrangements  of  the  ferrocinium  and 


, 
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picrate  ion; 

(2)  use  of  Cruickshank  and  MacDonald's^^  formula 
for  the  error  in  co-ordinates  due  to  the  wrong  hand 
when  using  anomalous  dispersion  indicated  that  both 
hands  were  incorrect  and  that  the  mean  -  a  disordered 
model  -  would  be  a  better  description;  and 

(3)  the  R  factor  was  consistently  higher  for  the 
non-centric  model  than  for  the  centric  model,  despite 
the  fact  that  the  latter  has  fewer  variable  parameters. 

Allowing  the  iron  atom  anisotropic  motion  and 
adding  the  hydrogens  of  the  cyclopentadienyl  rings 
as  hindered  rotors  with  C-H  distance  of  1.0A,  dropped 
the  R  factor  to  =  8.4%,  R2  =  12.3%.  Further  refine¬ 
ment  with  all  atoms  anisotropic  and  the  picrate  hydrogen 
added  gave  a  final  reliability  factor  of  R^  =  6.0%, 

R2  =  8.5%.  It  was  noticed  during  anisotropic  refine¬ 
ment  that  some  atoms,  noticeably  01  as  well  as  the 
disordered  nitro  group,  had  much  larger  thermal  para¬ 
meters  in  the  z  direction.  For  this  reason,  01  was 
also  moved  off  the  mirror  plane  and  given  a  half 
population  above  and  below  the  mirror  plane.  Refine¬ 
ment  of  this  disordered  model  gave  a  final  z  co-ordinate 
of  z  =  0.1794  ±  .1842.  Because  of  this  large  error  in 
the  z  co-ordinate,  01  was  returned  to  the  mirror  plane 
and  no  other  atoms  were  allowed  to  become  disordered. 
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In  the  final  cycle  of  refinement  no  shift  was  greater 
than  1/10  of  an  esd.  A  final  difference  map  gave  no 
peak  greater  than  0.42  e  /A?.  The  observed  and  calcu¬ 
lated  structure  factor  amplitudes  (x.10)  are  given  in 
Table  VIII.  Table  IX  gives  the  final  positional 
co-ordinates  of  the  individual  atoms  and  the  aniso¬ 
tropic  B  values  are  listed  in  Table  X.  Hindered  rotor 
parameters  for  the  cyclopentadienyl  rings  are  in 
Table  XI. 
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Table  VIII 


Observed  and  calculated  structure  factor  amplitudes 
in  electrons  for  [Fe (C5H5 ) 2+] [H2C6N307~] 
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0 

10 

302 

309 

0 

12 

188 

168 

1 

7 
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1 

11 
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2 

2 
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2 

4 
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299 
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10 

141 

207 

2 

12 
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3 

1 
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3 

7 
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4 

2 
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4 

4 
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5 

1 

216 
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5 

7 
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6 

4 
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6 

12 

168 

153 

7 

7 

184 

179 

7 

9 

193 

172 

L 

=  6 

0 

0 

392 

437 

0 

8 

172 

227 

1 

3 

188 

173 

2 

0 

262 

289 

2 

8 

161 

159 

4 

0 

178 

97 
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Table  IX 

Positional  co-ordinates  for  ferrocinium  picrate 


Atom 

X 

y 

z 

Fe 

0.0 

0.1851(1) 

0.2500 

Cl 

o 

• 

o 

0.4408(8) 

0.2500 

C2 

0.0932  (9) 

0.4829  (6) 

0.2500 

C3 

0.9464  (9) 

0.5500  (5) 

0.2500 

C4 

o 

• 

o 

0.5827  (7) 

0.2500 

Nl 

0.1996 (11) 

0.4511(5) 

0.2010(29) 

N2 

0.0 

0.6533  (7) 

0.2500 

01 

0.0 

0.3801(5) 

0.2500 

02 

0.2092(8) 

0.3960  (5) 

0.266(10) 

03 

0.2733  (9) 

0.4814 (5) 

0.1285(21) 

04 

0.0862(7) 

0.6822  (4) 

0.2500 

HI 

0.1680 

0.5730 

0.2500 

Parameters  without  estimated 


errors  were  not  refined. 


■ 

(V)i  •  • .  q 

Anisotropic  temperature  (B)  factors  for  ferrocinium  picrate 
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• 

(B)  fTC 


Table  XI 
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RESULTS 

A  drawing  of  the  ferrocinium  picrate  molecule  is 
shown  in  Figure  14.  The  packing  of  the  ions  projected 
onto  the  ab  plane  is  given  in  Figure  15.  The  atoms  in 
the  upper  picrate  ion  are  dotted  in,  while  those  below 
are  not.  Bond  angles  and  bond  lengths  are  given  in 
Table  XII,  while  some  interionic  contacts  are  listed  in 


Table  XIII. 


Figure  14 


Drawing  of  ferrocinium  picrate 
showing  numbering  system  used 


• 
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Figure  15 


Packing  of  ferrocinium  picrate 
as  projected  onto  the  ab  plane 


o 
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Table  XII 


(a)  Intramolecular  bond  lengths  in 

o 

ferrocinium  picrate  (A) 


Fe 

-  Cll 

2.074  (7) 

Fe 

-  C12 

2.070  (5) 

Fe 

-  C13 

2.064  (6) 

Cll 

-  C12 

1.395  (6) 

Cll 

-  Hll 

1.00 

H12 

-  01 

2.439  (8) 

H12 

-  02 

2.277  (10) 

01 

-  02 

2.641  (10) 

C2 

Cl 

• 

C2  ' 

(b)  Bond 

01 

— 

Cl. 

— 

C2 

Cl 

— 

C2 

— 

C3 

Cl 

— 

C2 

— 

Nl 

C2 

- 

C3 

- 

C4 

C3 

- 

C4 

— 

C3  ' 

C3 

- 

C4 

— 

N2 

C4 

— 

N2 

— 

04 

04 

- 

N2 

— 

04 

C2 

— 

Nl 

- 

02 

C2 

— 

Nl 

— 

03 

02 

— 

Nl 

— 

03 

dihedral  angle: 


Cl 

-  01 

1.231  (17) 

Cl 

-  C2 

1.445  (14) 

C2 

-  C3 

1.360  (15) 

C3 

-  C4 

1.357  (12) 

C2 

-  Nl 

1.517  (16) 

C4 

-  N2 

1.429  (18) 

Nl 

-  02 

1.210  (32) 

Nl 

-  03 

1.216  (15) 

N2 

-  04 

1.228(9) 

(C3 

-  HI) 

1.029  (11) 

angles  (deg) 

107.6  (14) 

126.2  (7) 

126.9  (12) 

117.2  (9) 
118.5(12) 

121.5  (14) 
119.2(7) 

118.6  (8) 

122.9  (15) 
113.  (2) 
123.  (1) 
123.  (2) 

(z  =  0.25  plane)  and 


(Nl, 


02,  03  plane)  26.6° 


• 

■ 

' 

- 
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Table  XIII 
Interionic  contacts 


Fe 

-  01 

3.951(11) 

03 

-  H13 

2.383  (12) 

01 

-  Hll 

2.439  (18) 

03 

-  H13 

2.668  (12) 

01 

-  02 

2.641(10) 

03 

-  03 

1.928  (29) 

01 

-  Cll 

3.234  (12) 

04 

-  Hll 

2.505(7) 

01 

-  04 

3.830(6) 

02 

-  Hll 

2.277  (10) 

■ 
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DISCUSSION 

The  compound  ferrocinium  picrate  is  an  ionic  species 
consisting  of  the  ferrocinium  ion  Fe(Cp)2+  in  which  the 
iron  may  be  considered  to  be  in  the  formal  +3  oxidation 
state,  and  the  picrate  ion  [H-CrN-0_ ]  " . 

The  ferrocinium  ion  has  two  cyclopentadienyl  rings 
which  are  related  by  a  mirror  plane,  and  thus  unlike 
ferrocene  in  the  solid  state35'36,  the  rings  are  eclipsed. 
Ferrocene  in  the  gaseous  state  6,  as  well  as  ruthenocene33 
in  the  solid  state,  are  also  eclipsed.  The  Fe-C  distances 

o 

were  found  to  be  longer  (average  2.070 (5) A)  than  those  of 
36  ° 

ferrocene  (2.047 (5)A)  while  the  average  C-C  distance  of 

°  o 

1.395 (6)A  for  ferrocinium  was  shorter  (1.416 (4)A)  as 
obtained  from  the  hindered  rotor  model. 

Although  distortions  from  symmetry  may  be  present, 

because  of  the  disorder  in  the  crystal  and  large  oscilla¬ 
tions  of  the  rings  (RMS~18°),  the  results  of  any  indi¬ 
vidual  atom  refinement  would  likely  be  unreliable.  Evi¬ 
dence  of  this  can  be  seen  in  the  results  of  Pettersen^  in 

o  o 

which  the  C-C  distance  ranges  from  1.33A  to  1.44A  with  an 

o 

error  of  +0.02A  or  greater.  Use  of  the  hindered  rotor 
model  allows  for  the  oscillation  of  the  rings  and  should 
give  a  good  average  value  for  the  C-C  distance.  The 
lengthening  of  the  Fe-C  distances  and  shortening  of  the 


1 

' 

■ 

m 
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C-C  distances  in  the  ferrocinium  ion  is  in  agreement 

4  2  43 

with  the  molecular  orbital  calculations  '  and  spectral 
44-48 

results  which  indicated  that  the  ground  state  of 

ferrocene  in  the  eclipsed  formation  is  E2g  ^aig^  2  ^e2g^ ^ 


and  that  of  the  ferrocinium  ion  is 


E2g[(alg)2(e2g)3]- 


Since  a^  and  e ^  are  the  highest  bonding  orbitals, 
removal  of  one  electron  from  the  e~  level 

2g 


e . 


ig 


e .  * 
-  ig 


-H— aig 

ferrocene 


•H-H 519 

ferrocinium  ion 


would  weaken  the  Fe-C(ring)  bonds  and  result  in  a  longer 
Fe-C  distance  in  ferrocinium  ion  than  in  ferrocene.  It 
would  also  strengthen  the  C-C  bond  thus  shortening  it  in 


' 
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FeCCp^*  since  an  electron  is  being  removed  from  a  C-C 

42 

antibonding  orbital.  These  calculations  were  done 

assuming  symmetry.  However,  the  symmetry  of  the 

ferrocinium  ion  in  this  structure  is  D_,  rather  than 

5h 

.  Molecular  orbital  calculations  using  D,_^  symmetry 
give  the  same  pattern  of  energy  levels.  If  the  assump¬ 
tion  is  made  that  the  order  of  the  energy  levels  is  the 
same  in  the  eclipsed  form  as  in  the  staggered  form, 

then  a  ground  state  of  the  ferrocinium  ion  would  be 
2  3 

E2'  ^al'  '  ^2  ^  an<^  t*ie  ar9ument  given  above  would 

still  apply. 

The  picrate  ion  lies  on  a  mirror  plane  at  z  =  0.2500 

and  across  a  second  mirror  at  x  =  0.0.  However,  three 

atoms,  namely  those  of  the  ortho  nitro  group  (Nl,  02 

and  03),  were  allowed  to  move  off  the  mirror,  and  have 

"half  atoms"  above  and  below  the  z  =  h  mirror.  The 

reasons  for  this  were  discussed  in  some  detail  earlier. 

Although  because  of  the  disorder  problem  the  results 

are  somewhat  suspect,  the  geometry  calculated  from  the 

disordered  picrate  model  agrees  fairly  well  with  those 

39 

of  potassium  and  ammonium  picrate  .  The  N-0  distances 

o  o 

within  the  ortho  nitro  group  (1.210 (32) A  and  1.216  (15) A 
are  the  same  within  experimental  error  as  those  of  the 

o  o 

potassium  and  ammonium  structures  (1.232 (5) A,  1.229 (5) A, 

o  o 

1.237 (5) A  and  1.206 (5) A)  respectively.  The  ortho  nitro 
group  was  found  to  be  twisted  26.6  degrees  from  the 


. 
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plane  of  the  benzene  ring,  while  in  the  potassium  and 
ammonium  compounds  it  was  25°  -  27°.  The  C-N  bond 

o 

is  slightly  shorter  (1.429 (18) A)  than  those  found  by 

o  o 

Maartmann-Moe  (1.457 (6) A  and  1.461 (5) A).  The  rest 
of  the  molecule  also  has  a  geometry  very  similar  to 
the  potassium  and  ammonium  picrates. 


Chapter  III 


THE  CRYSTAL  STRUCTURE  OF  HYDRIDOTRIPHENYLSILYL 
(tt-CYCLOPENTADIENYL)  DICARBONYL  MANGANESE 


79 


INTRODUCTION 

The  ultraviolet  irradiation  of  metal  carbonyls 

in  the  presence  of  substituted  silanes  has  lead  to 

the  preparation  of  a  new  class  of  metal-silicon  hydrides.. 

These  compounds  seem  to  fall  into  two  groups;  for 

example,  those  prepared  from  triaryl  silanes  such  as 

(tt-C^H,.)  (OC)  2Mn  (H)  Si  (C^H^)  (1)  ,  and  those  containing 

the  more  electronegatively  substituted  silyls  such  as 

(tt-C^H^)  (OC)  2Mn  (H)  SiCl^  (2).  The  latter  group  of 

compounds  shows  well  defined  infrared  bands  in  the 

region  of  "normal"  terminal  metal-hydrogen  stretching 

frequencies  (1887  cm  \  deuteride  1360  cm  The 

former  group,  however,  showed  no  bands  in  either  terminal 

4 

or  bridging  M-H  regions.  The  solid  state  Raman  spectra 
shows  a  weak  broad  band  at  1900  cm  ^  which  shifts  to 
1350  cm  ^  in  the  deuteride,  these  bands  being  in  what 
is  considered  the  terminal  metal-hydrogen  region.  In 
order  to  understand  the  differences  in  the  spectral 
results,  the  crystal  structures  of  (1)  and  (2)  were 


undertaken. 


. 


. 
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EXPERIMENTAL 

The  compound  hydridotrichlorosilyl  (iT-cyclo- 

pentadienyl)  dicarbonyl  manganese,  (n-CirHt- )  (OC)  _ 

5  5  2 

Mn(H)SiCl^  (2),  forms  light  brown  needle  crystals. 
Rotation,  VJeissenberg  (h0i-h4i  and  0k£-9k£)  and 
precession  (h0£-h3£)  photographs  using  several 
crystals  showed  the  crystals  to  be  monoclinic  with 
unit  cell  a  =  14.71A,  b  =  11.46A,  c  =  14.85A  and 
3  =  107.75°.  Systematic  absences  of  hk£;  h+i  = 

2n+l ,  OkO;  k  =  2n+l  and  h0£;  &  =  2n+l  fixed  the  space 
group  at  B21/c  (a  non-standard  setting  of  P21/c) . 
Eight  molecules  per  unit  cell  gave  a  calculated 
density  of  1.75  gm/cc.  Since  the  molecule  does 
not  possess  a  centre  of  symmetry,  it  is  most  likely 
that  the  molecules  lie  in  general  positions  with  one 
molecule  per  asymmetric  unit.  The  crystal  was 
mounted  on  the  PAILRED  automatic  diffractometer  in 
order  to  collect  data.  However,  measurement  of 
standard  reflections  showed  that  decomposition  of 
the  crystal  in  the  X-ray  beam  was  so  rapid  (over 
30%  in  20  hours) ,  an  accurate  crystal  structure  was 
impossible.  The  determination  was  terminated  at 
this  point. 


'  . 
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The  compound  hydridotr iphenylsilyl  (Tr-cyclo- 

pentadienyl)  dicarbonyl  manganese  (7T-CrHr)  (OC) 

5  5  2 

Mn (H) Si (C6H5)  was  prepared  by  J.K.  Hoyano49  by 
the  ultraviolet  radiation  of  (C5H5)Mn (CO)  in  the 
presence  of  an  excess  of  triphenylsilane  and  re¬ 
crystallized  from  a  mixed  hexane-dichloromethane 
solvent  system  to  obtain  crystals  suitable  for  single 
crystal  X-ray  diffraction  study.  The  yellow  crystals 
were  mounted  parallel  to  an  extinction  direction  (a 
face  diagonal) .  Rotation,  Weissenberg  (0k£-2k£) 
and  precession  (kh0~hk2  and  h0£-h2£)  photographs 
showed  the  crystal  to  be  monoclinic.  Systematic 
absences  of  OkO,  k  =  2n+l  and  h0£,  h+ £  =  2n+l 
fixed  the  space  group  as  P21/n.  The  general 
positions  in  this  space  group  are  x  y  z 

x  y  z 
^+x,  h~y f  h+z 
h-x,  h+y ,  h-z 

The  unit  cell  dimensions  were  obtained  from  least 
squares  refinement  of  20  values  of  reflections  that 
had  been  carefully  centred  in  the  counter  window  on 
the  Picker  four  circle  manual  diffractometer,  using 

o 

CuKap  radiation  (A  =  1.5405A).  The  cell  dimensions 
are  a  =  13.200(3)A,  b  =  17.328(4)A,  c  =  9.438(2)A 
and  3  =  92.47(4)°.  The  density  was  determined  by 


' 
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flotation  to  be  1.347(3)  gm/cc  which  agrees  well  with 

the  calculated  density  of  1.342  gm/cc  (Z  =  4, 

°  3 

V  =  2156. 8A  )  and  is  consistent  with  one  molecule  per 
asymmetric  unit. 

The  crystal  used  in  data  collection  was  a  prism 
bound  by  the  faces  010,  010,  203,  203,  212  and  2l2 
with  dimensions  0.17  mm  x  .097  mm  x  0.097  mm.  It  was 
mounted  along  a  face  diagonal,  the  a*  axis  being  co¬ 
incident  with  the  diffractometer  0  axis.  Data  was 
collected  manually  on  the  Picker  four  circle  diffracto¬ 
meter  using  CuK  X-radiation  and  a  coupled  w/20  scan 
technique.  The  angular  settings  (0,  x*  29)  were 
calculated  using  MIXG2  for  unique  sections  of  the 
reciprocal  lattice  limited  by  sin  0/A  ,<0.5.  A 

scan  width  of  2  degrees  at  a  scan  speed  of  2  degrees 
per  minute  was  employed  with  a  stationary  background 
count  of  30  seconds  on  each  side  of  the  scan.  The  beam 
was  monochromated  by  means  of  a  graphite  crystal  (002 
reflection) .  The  diffracted  beam  collimator  aperture 
was  0.5  mm  and  the  pulse  height  analyser  was  centred 

on  the  CuK  peak.  Intensities  were  calculated  using 
a 

I  =  P-t(B1+B2)  where  P  is  the  peak  count,  B1  and  B2 
are  the  background  counts  and  t  is  the  ratio  of  total 
time  of  peak  count  to  total  time  of  background  count. 
Two  thousand,  two  hundred  and  thirteen  unique  (hk£ 


■ 


. 


. 
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and  hk£ ,  h,  k  and  £  positive)  reflections  were  measured 

of  which  662  were  rejected  using  the  criteria  (1)  1^0, 

or  (2)  I^3a  where  a  =  [P  +  t^ (B^  +  B^)]  2.  Lorentz  and 

polarization  corrections  were  applied  to  the  remaining 

1551  reflections.  Twelve  well  distributed  standard 

reflections  measured  at  various  times  throughout  the 

data  collection  indicated  no  decomposition  of  the 

crystal.  The  maximum  variation  was  ±5  sigma.  There 

was  no  apparent  pattern  in  the  peak  height  variation. 

Absorption  corrections  (y  =58.8  cm-1)  were  applied 

c  u 

to  the  data.  Transmission  factors  ranged  from  0.51 
to  0.63  with  the  majority  lying  close  to  0.59. 


. 
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SOLUTION  AND  REFINEMENT 

The  three  dimensional  Patterson  map  gave  a  set 
of  three  peaks  agreeing  with  a  manganese  position  of 
x  =  -0.112,  y  =  0.250  and  z  =  0.060.  A  further  set 
of  peaks  gave  the  position  of  the  silicon  atom  as 
x  =  -.043,  y  =  0.199  and  z  =  -.165.  A  structure 
factor  calculation  using  these  co-ordinates  returned 

R1  =  E 1 1 Fo | - [ Fc 1 1  of  51.6%  and  R2  = 

I  |Fo  | 

of  64.0%.  The  weights  used  were  calculated  from 

a  (F)  s  w  =  1/a  (F)  2  where  a  (F)  =  [P  +  t2  (B  +B2)  + 

2  ^ 

p  I ]  2 ,  p  =  .03.  The  scattering  factors  were  calculated 

50 

from  Cromer's  co-ef ficients  for  all  atoms  except 

51 

hydrogen.  These  were  obtained  from  Mason  and  Robertson 
A  difference  map  phased  using  the  manganese  and  silicon 
atom  positions  revealed  the  location  of  the  two  carbonyl 
groups.  A  further  difference  map,  phased  using  the  six 
atoms  located  gave  the  positions  of  all  remaining  non¬ 
hydrogen  atoms.  After  one  cycle  in  which  the  phenyl 
groups  were  treated  as  rigid  bodies  and  the  cyclopenta- 
dienyl  ring  as  a  hindered  rotor,  the  R  factor  dropped  to 
R^  =  28.6%  and  R2  =  38.9%.  The  starting  parameters  for 
the  four  rings  were  calculated  using  MMMR  with  a  radius 

o 

of  1.397A  for  the  rigid  bodies.  At  the  end  of  iso¬ 
tropic  refinement  the  R  factor  was  R^  = 


Sw 1 | Fo | - | Fc | 

i  .2 

Zw  Fo 


8.6%  and 


. 


. 

.  >li  "t 
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1*2  =  10.5%.  A  difference  map  at  this  time  showed 

the  hydrogen  atoms  attached  to  the  aromatic  rings, 

.  —  °  8  —  °  8 
the  peak  heights  ranging  from  0.76e/A  to  0.38e/A  . 

These  were  added  as  rigid  bodies  with  the  same 

orientation  and  centre  of  gravity  as  the  carbon  rings 

to  which  they  were  bonded  and  with  a  C-H  distance  of 

o 

1.0A.  The  hydrogens  on  the  cyclopentadienyl  ring 

were  added  as  a  hindered  rotor  rather  than  as  a  rigid 

body.  The  difference  map  also  showed  a  peak  in  a 

bridging  position  between  the  manganese  and  silicon 

atoms.  This  peak  was  of  the  same  magnitude  as  those 

—  °  3 

of  the  twenty  aromatic  hydrogens  ('-0.5  e/A  )  and 

was  the  twelfth  largest  peak  on  the  map.  Adding 

the  hydrogen  atoms,  giving  the  manganese,  silicon 

and  carbonyls  anisotropic  temperature  factors  of 

2  2  2 

the  form:  exp-fp^h  +  322k  +  ^33^  +  $12hk  + 

$23kj^  and  allowin<?  for  anomalous  disperson 

(Af,MN  =  "-050'  Af"Mn  =  2’90'  Af’si  =  °-20  and 
Af"  .  =  0.40)  dropped  the  R  factor  to  R  =  6.2%, 

S  X  ± 

R2  =  7.3%  after  four  cycles.  A  difference  map  showed 
that  the  peak  previously  found  at  the  bridging  posi¬ 
tion  (x  =  -.009,  y  =  0.248  and  z  =  -.014)  was  still 

—  0  3 

present  with  a  peak  height  of  0.53e/A  and  was  the 
second  largest  peak  present.  This  hydrogen  peak  is 
of  very  similar  density  to  those  of  the  hydrogen  peaks 
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found  by  Ibers  et  al .  in  the  structure 
CoH(N2)  (P  (CgH,. )  ^)  3  .  The  largest  peak  was  located 
near  the  para  carbon  of  a  phenyl  ring  (C44) ,  and 
since  it  disappeared  on  the  low  angle  (sin  0/1^0.30 
difference  map,  was  most  likely  due  to  the  thermal 
motion  of  the  phenyl  ring.  To  further  check  the 
validity  of  the  hydride  peak,  difference  maps  were 
calculated  using  various  sin  6/A  data  cutoffs 
(sin  e/A^0.35,  ,$0.30,  ^0.25,  and  .$0.20).  In  all 
four  maps,  the  hydrogen  peak  was  the  largest  peak  on 
the  map  and  was  approximately  a  factor  of  two  greater 
than  the  next  highest.  The  peak  positions  remained 
constant  on  all  four  maps  while  its  height  varied 
according  to  the  sin  0/A  cutoff.  Table  XIV  gives 
the  observed  and  calculated  peak  heights.  The 
calculations  are  the  same  as  those  of  La  Placa 
and  Ibers^ . 

The  hydridic  hydrogen  was  added  as  an  isotropic 

atom  and  allowed  to  refine.  At  the  same  time  the 

carbons  in  the  phenyl  rings  were  refined  using  the 

individual  atom,  anisotropic  temperature  factor  model. 

Three  cycles  of  refinement  gave  a  final  R  factor  of 

R  =  4.2%  and  R2  =  4.7%.  In  the  final  cycle  no 

shift  was  greater  than  1/10  of  an  e.s.d.  A  final 

_  0  3 

difference  map  showed  no  peak  greater  than  0.19e/A 
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Table  XIV 

Hydrogen  peak  as  a  function  of  sin  0/A 


sin  6/A 
cutoff (A“l) 

Number  of  terms 
in  the  unique  section 

F1  __  °  o 

P  (e/AJ) 

c  ,  -7°  3  x  b 

P  (e/A  ) 

0.20 

149 

0.20 

0.16 

0.25 

278 

0.29 

0.25 

0.30 

464 

0.37 

0.34 

0.35 

703 

0.41 

0.42 

0 . 50a 

1551 

0.53 

0.53 

a  full  data  set 


c  r  s 

p  =  1 


2tt 


o  ,,  .  2  2,., -2  ,  ^  2  ,,  r  2.  2  , 

(1+a  s  /4)  exp(-Bs  /16tt  )s  ds 


s  =  4nsin  0/A 
a  =  Bohr's  radius 
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or  less  than  -0.20e/A  . 

Table  XV  lists  the  observed  and  calculated 
structure  factor  amplitudes;  10|Fq|  and  10|F  |, 
both  in  absolute  units  of  electrons.  The  final 
positional  parameters  for  the  individual  atoms 
are  given  in  Table  XVI  and  their  anisotropic  temp¬ 
erature  factors  in  Table  XVII.  Parameters  for  the 
hindered  rotors  are  found  in  Table  XVIII.  The 
estimated  standard  deviations  were  obtained  from 
the  inverse  matrix  of  the  final  least  squares  cycle. 


. 

. 


- 


89 


Table  XV 

i 

Observed  and  calculated  structure  factor  amplitudes 

(x  10)  in  electrons  for  (tt-C^H  )  (OC)  0Mn  (H)  Si  (C^H_ )  _ 

->  ->  l  6  5  3 


0  576  5*J 

r  4?s  4i6 

9  797  ?38 

o  40*  tr» 

I  1/?  1 1 T 

I  108?  110? 


n  y 
?  oj 

700 

176 

399 

*09 

410 

445 

750 


197 

730 

658 


43* 

707 

784 

837 

437 

716 

115 

4«6 


176 
346 
359 
179 
359 
48? 
1 150 
70* 

no 

748 


793 

400 

487 

75C 


770 

665 

547 


838 

4*7 
704 
119 
4  4* 
75? 
757 
108 


157 

398 

47* 


610 

585 

714 

414 


616 

594 

77* 

94« 

171 

767 

310 

754 


488  489 


779 

719 

759 


473 

775 

774 

355 

10" 


718 

418 

763 

706 

747 

71 

303 

744 

*36 

316 

74? 

7*3 


1*0 

135 

774 

*73 

770 


77? 

7C4 

137 


739 

130 

175 


1103 

3*1 

90* 

657 


1791 

351 

943 

657 

137 

674 

70*7 

*73 

361 


984 

1071 

900 

1569 


139 
779 
609 
I  101 

*99 

1848 


157 

318 

39* 

754 

784 


187 

264 

914 


7*0 

292 

398 

979 


1125 

7*9 

294 

4*2 

•o? 


1237  *  II 

2*5  -1  1| 

7*3  -7  1 1 

*7?  -3  11 

195  -4  11 


?»■> 

703 

600 

515 


777 
2  42 
139? 


130 

798 

778 

1*11 

2C7 


104 

•5 

987 

71? 

740 

15? 

094 

130 

6*3 

307 

7C0 

1017 

104 

162 

340 


377 

735 

386 

269 


559 

607 

644 


1*9 

*90 

l*o 

310 

961 

*65 

713 

?** 

6«4 

386 

19* 

71* 

746 


627 
465 
32* 
*7? 
1*7 
174 
170 
4*0 
98 
*34 
394 
54  8 
160 
343 


306 

15? 

415 


270 

943 

717 

113 

7*9 

*45 

*73 

745 

*04 

537 


700 

31? 

710 

1039 


640 

455 

114 


423 

677 

346 

223 

366 


30* 

444 

310 

744 


431 

511 

64* 

255 

130 

476 

251 

*76 

77* 


777 

371 

786 

1*7 

17* 

19« 


564 

211 

386 

2*6 


253 
202  ' 
244 
16* 


202 

306 

700 

3*0 


287 

777 

964 


179 

*18 

476 

3*5 

14? 

108 

230 

217 

270 

129 

71* 

103 


*19 

*34 

378 


257 

718 

*05 


209 

504 

476 

95 

44* 

708 

716 


108 

3*0 

*16 


*78  468 

601  59* 

770  77* 

567  534 

770  764 

2*8  294 

1049  1103 

*83  48? 

294  316 


411 

71) 


6*3 

823 

39* 


825 

705 


1127 

644 

1*0 


107 

414 

137 


3*3 

284 

*37 


312 

670 

276 

*90 

375 

270 

630 

2C0 

35* 


186 

7*7 

314 

565 

630 

360 


485 

497 

7*8 

106 

2*7 

31* 

356 

305 

717 

766 

468 

1*4 

199 

111 

412 

656 

6)8 

221 

340 

444 


447 

*77 

668 

22  7 


768 

547 

7*6 

*7* 


710  143 

730  ?»| 

1295  1747 

344  340 

741  291 

III  95 
975  471 

270  759 


176  129 

76*  27) 

401  4*6 

192  1*6 

196 
18* 

775  ?70 

*0?  599 

2*9  234 

7)6  24? 

2)0  225 


168 

181 


187 

152 


1*1 

194 


274  2?  T 

4*7  *57 

254  263 


7*6  234 

3*3  3*7 

36*  34  7 

149  1*6 

21)  719 

189  187 

1*5  151 

165  172 

179  |61 

7*  84 

1014  104? 


107 

*71 

304 

419 

717 


113 

*31 


14* 

115 
654 
315 
*24 
749 
164 
76*  257 

715  272 

221  705 

19*  2C5 

765  766 

4)6  4)3 

508  *7? 

99  113 

242  254 

138 
645 
378 
348 
1*6 

716  ?0t 

7*8  255 

177  114 


225  206 

9*  105 

14)  136 

170  1*0 

220  217 

765  29) 


705  214 
223  219 
266  768 


141 


136 


335  137 

560  458 

49?  444 

217  213 

197  146 

364  344 

22 *  211 
14J  137 

96  56 

170  170 


2  12  423 


-l  12  352 


98 


189 

*7 

444  470 

38?  174 

243  306 

164  172 

365  36) 

345  353 

158  152 

386  380 

123  120 

217  287 

267  273 

3*5  401 

106  6* 
275  261 

116  109 


170 


141 

155 


111 


111 


1)2  118 
89*  *44 

20*  207 

765  7*8 

1641  1762 

612  626 
724  752 

178  It* 
1)3  112 

217  213 

717  7*0 

5*4  548 

312  308 


72 


50 


144 

124  372 
354  355 
1*4  175 
514  513 


9*7 


335 


1*7 

8* 


624  08 

*71  *74 

827  891 

38*  386 

182  |9) 

127 
41 
80 

161  173 

*20  407 

*84  4  75 

24?  791 

208  207 

3*8 
17? 


187 

334 


244 

109 

484 

177 


228 
355 
275  290 


164 

403 


»7| 
147 
190 
149 
211  218 
124  132 

228  220 
385  150 

186  176 

*16  8*4 

344  383 

147  123 

29*  104 

240  238 

38?  357 

57?  572 

377  378 

274  783 

360  35? 

173  37* 

130  143 

414  413 

3*4  3)6 

170  174 

It)  184 
762  264 

557  575 

484  47) 

3)1  336 

270  260 

*54  *57 

2  77  277 

21*  214 

124  123 

103  107 

48)  455 

*37  446 

42 


71 

356 


350 
160 
274  777 

22)  214 

259  ?5l 

140  154 

778  781 

411  402 

416  616 

33*  374 

437  446 

244  274 


100 


84 


754 


314 
757 
197  198 

150  154 

474  664 

514  520 

241  232 

115  107 

4)3  449 

180  182 
235  214 

2*9  253 

1)4  190 

84  63 

3  0*  29* 

124  135 

340  328 


926 
2)8 
144 
121  111 


152 

357 


1)8 

40 


153 
355 
270 
177  182 

468  444 

261  269 

37*  403 

24*  247 

42  41 

1)5 
118 
433  447 

518  906 

554  543 

145  1*9 

7*0  244 

466  494 

7*0  771 

154  147 

247 
10* 
1*4 
1 7?  1*3 

1)6  140 

4 14  406 

220  214 

34?  93* 

240  28? 

429  438 

142  172 

152  158 

177  173 

146  148 

100  41 

154  148 

441  4*9 

84  4* 

371  37? 


120 

174 


36* 

170 


381 

144 


194 
3*0  371 

627  644 

10)  110 
414  414 

912  91) 

2  T*  294 
22? 


41 

184 

154 


47 

700 

144 


164  148 

21?  211 
1*4  |84 

264  280 

597  597 

lit  925 
122  125 

414  494 

3*4  980 

46*  4*2 

1*5  |4S 

22*  220 
119  *6 

2*4  299 

2  7*  2*7 

*5*  *52 

20*  207 

3*1 


141 


134 
43  44 

m  io7 
115  127 

111  30* 


147 

104 

194 

190 


369  334 
147  1*6 
226  713 
276  217 
1)4  150 
136  160 
240  220 


144 

214 


161 

71* 


7)4 

84 

?0« 


7*4 
47 
?C4 
347  367 

5*0  576 

459  459 

IK4  1141 
134 
760 
7*7 
27)  26* 

607  598 

20?  209 

376  310 

15*  l*o 
149  |86 

7*  84 

160*  1?6 


176 

7*7 

7*4 


I  93 


777 

II* 


*84 


40 

774 


102 


54 


116 

14? 


247  258 
447  48* 
674  *09 


178 


II® 


3*4  960 
37®  3?o 
770  718 
210  710 

1*9 


103 


84 


447  424 
248  748 
767  266 
714  219 
674  596 
479  559 


’99 


398 


787  263 
2®5  746 
133  114 


154 

160 


148 

164 


?*• 


137 
74* 
701  704 

*53  46) 

75"  273 

477  43? 

174 


393 


961 


111  171 

122  125 

40  66 

*78  43® 

171  171 

*7)  *70 

754  ?42 

506  *88 

40?  )•• 

14?  704 

337  3)4 

217  ?48 

119 


112 

307 


764 

456 

110 

145 

13* 

1*9 

463 


1)4 

176 

36* 

16) 

300 

346 

710 

70) 


150 
■  7  80 

?"*  774 

714  715 

20  2  715 

2»»  246 

T|»  73? 

?)*  250 

1°) 


771 
1  35 


70* 

618 

387 


65  14 

736  778 

*06  506 

748  77® 

13®  144 

?">  786 

*03  497 


14? 


10? 

170 
517 
49*  544 

774  260 

515  515 

’47  394 

*1  44 

164  16* 


470 


1"5  |I8 

14?  211 

75®  7*9 

116  ID 
404  *10 

677  434 


89 

770 


233  778 


10® 
450  454 

314  )05 

346  340 

*15  412 

275  763 

161 


154 

?3C 

109 

1*4 

2M 


440 

1»* 

??4 

790 

176 

Ml 

12® 


563 

733 

117 

146 


314  304 
541  59? 
476  579 
>69  ?"4 

331 


215 


770 
178 
4*9 
7® 1  289 

?6l  766 

4’S  409 

34’  347 


765 


747 


768  766 
279  725 
754  234 


53*  444 

_  1 39  1*1 

100  5? 

156  139 

*®6  *09 

*1  3  394 

573  *97 


84 


43 


314  317 


431 


*31 


19* 

737  256 

24*  297 

639  6*0 

31  3 


341 
6*  7 


377 

665 


213  7C7 

7*3  744 

354  374 

237  735 


171 


174 
®3  75 

777  214 

*75  415 

315  309 

238  744 

27?  28C 

7)4  757 

35«  3*1 

*98 
398 
?4?  276 

'6?  369 

2")  146 

119  1?4 

*47  47? 

412  41) 

164 
11« 

727 
30? 

264  ?56 

707 
1)8 


378 


168 

119 

23? 


387 

90 


7? 

360 


3  "  1 
9* 

147 
1*7 
13« 

84 
41 
372 
167 
385  ’64 

100  181 
127  88 

167  |96 

149  145 

761  270 

30?  744 

453  46) 

527  518 

14)  I 4 | 

275  278 

164  168 

347  937 

170  104 

?0?  196 

771  76) 

1)1 
264 
703 
?*?  755 


110 

796 


194  177 
214  216 
7)9  779 
2)9  230 


-1  12  174 


111 

1*8 


041  424 

140  203 

2*0  26) 
114  90 

II?  44 

143 


147 

2*2 


236 

12? 

15* 

7)7 


124 

354 

II* 


127 

167 

70S 

166 

149 

68 

1<*7 

19* 

156 

127 

216 

70? 

14) 

144 

143 

175 


3*0  347 
76"  764 
46?  *67 
?®6  )1| 
777  714 
169  16) 
491  450 
356  372 
448  460 
275  277 
112  IC3 
45)  458 
130  123 


303  906 

120  103 

504  *8* 

234  2)3 

517  429 

1)0 


335 


34? 


140  171 

134  107 

254  25* 


47 


77 


117  123 

IT)  197 
106  122 
*9)  53* 

253  768 

178  1)6 

172  168 

980  366 


128 


12? 


349  3)6 

II*  1)0 
271  2*9 

244  7)4 

1*5  206 

157  157 

194  II? 
147  145 

215  20* 

157  149 

184  1*3 

17*  115 


191 


17* 


142 

288 


150 


1*7 


166 

457 

10* 


70* 

1*1 


216 
1®* 

136  171 

*62  *70 

267  7*3 

204  206 

705  202 

7C3  183  ' 

75?  752 

m  i*7 
172  II* 
l«l  195 
44*  4*7 

270  765 

380  39* 

1)3  17* 

917  31* 

77)  776 

751  2*5 

142  190 

3*1  37) 

171  90 

10*  *9 

23*  275 

1*7  1)7 

98  120 

71*  215 

1)9  IS* 
372  33* 

77)  764 

14*  1*7 

21*  772 

Jfc)  37* 

28 4  300 

)*>  3?) 

778  773 

22*  234 

4*6  *40 

77*  27? 


18*  20) 
162  143 

304  313 

297  29* 

1)7  13* 

11*  17* 

299  311 

151 
294 
1)9 

17* 

199 
177 
718  717 

*76  47? 

177  171 

168  166 
79*  796 

1*1  197 

117  11? 

7)7  2)7 

170  173 

49  94 

371  304 

10)  97 

13*  13* 

127  10) 


161 


170 


363  3)3 

101  102 
307  31) 

210  21* 
1*0  17* 

134  149 

1C8  107 


179 


171 


269  25* 

125  '  9* 

171  173 

2*1  246 

217  278 

107  111 

7*3  2)9 

149  14* 

10*  100 
742  238 

19*  1*| 

12*  117 

125  108 

1*5  17* 

107  94 

340  334 

17)  157 

400  40) 

189  167 

1)1  137 

20)  200 
144  160 

320  317 

1)6  119 


103 


91 


158  130 

29*  291 

370  920 

197  197 

?43  269 

151  173 

277  763 

15)  157 

170 
70 
1)8 
1 37 
7)5 


91 

146 


200 

271 


21* 

27* 

770 

297 

290 

2)9 


216 

1C* 


2)9 
9*5 
790 
713 
127 
1C* 
710 
761 
739 
78* 
777 
7)4 
176 
I  IT 
145 
182 

153 
165 
396 

154 
127 
71* 


160 

375 

112 

217 


16* 
37*  31* 

19)  It* 
*0 1  60* 

1*1  196 

4?0  4)7 

717  712 

149  16* 

200  14# 

2  >0  279 

2 CO  179 
1*5  1)7 


291 


21« 

160 

248  7*9 

77*  766 

341  3*0 

770  7*7 

11*  173 

728  736 

1*8  156 

116  l«l 
719  29| 


307 
17* 
3)0 
2  35 
2)6 
1)6 
108 
II* 
257 


153 

1*3 

215 

17* 

529 

160 

II* 

6*2 


))• 

154 

267 


31* 
30* 
126 
33* 
72* 
.  23* 
123 


228 

1*3 

37* 


2*1 

22* 

207 

14) 


1*9 

43* 

120 


11* 

178 

101 

30* 

II* 


167 

288 

22) 


13* 

I2T 

708 

13* 

109 

37* 


187 

189 


179 


■ 
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Table  XVI 

Final  atom  co-ordinates  for  (CcH_ )  (OC) 0Mn (H) Si (C.tU ) _ 

b  b  Z  obi 


Atom 

X 

y 

z 

Mn 

-0.1187  (1) 

0.2382(1) 

0.0658  (1) 

Si 

-0.0492(1) 

0.1991(1) 

-0.1567(1) 

Cl 

-0.0827(4) 

0.3331(4) 

0.1165 (6) 

01 

-0.0578(3) 

0.3956  (3) 

0.1474  (5) 

C2 

-0.2109(4) 

0.2782(3) 

-0.0535(6) 

02 

-0.2751(3) 

0.3056(2) 

-0.1262(5) 

C21 

-0.1459 (4) 

0.1368  (3) 

-0.2546  (6) 

C22 

-0.2167(5) 

0.1689  (3) 

-0.3513 (6) 

C23 

-0.2910(5) 

0.1250  (5) 

-0.4210  (7) 

C24 

-0.2964  (6) 

0.0476  (5) 

-0.3957(8) 

C25 

-0.2278(7) 

0.0136(4) 

-0.3015(8) 

C26 

-0.1538  (5) 

0.0574(4) 

-0.2314 (6) 

C31 

-0.0206(4) 

0.2819(3) 

-0.2787(6) 

C32 

-0.0198  (4) 

0.3590  (3) 

-0.2374  (6) 

C33 

0.0085  (5) 

0.4173(3) 

-0.3274  (7) 

C34 

0.0374(4) 

0.4005  (3) 

-0.4615(7) 

C35 

0.0382  (4) 

0.3247  (4) 

-0.5066  (6) 

C36 

0.0096(4) 

0.2661(3) 

-0.4157 (6) 

C41 

0.0721(4) 

0.1415  (3) 

-0.1424  (6) 

C42 

0.0997(5) 

0.9020  (3) 

-0.2512  (6) 

C43 

0.1950(6) 

0.0545(4 

-0.2480(8) 

C44 

0.2608(5) 

0.0676(4) 

-0.1333 (10) 

C45 

0.2378  (5) 

0.1154(4) 

-0.0260(8) 

C46 

0.1435(4) 

0.1515  (3) 

-0.0306(7) 

HI 

-0.0132(32) 

0.2467(23) 

0.0020  (45) 

• 

Anisotropic  temperature  factors  (x  10 
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Table  XVIII 
(Contd. ) 


(b)  Derived  positional  parameters  for  hindered 

rotors  and  rigid  bodies 


Ring 

Atom 

X 

y 

z 

1 

H22 

-.2128 

.  2274 

-.3701 

H23 

-.3415 

.1510 

-.4899 

H24 

-.3505 

.0152 

-.4451 

H25 

-.2308 

-.0440 

-.2805 

H26 

-.1021 

.  0324 

-.1607 

2 

H32 

-.0427 

.  3657 

-.1358 

H33 

.0051 

.  4708 

-.2833 

H34 

.0566 

.4465 

-.5192 

H35 

.0603 

.3171 

-.6076 

H36 

.0126 

.2119 

-.4600 

3 

H42 

.0490 

.0831 

-.3347 

H43 

.2109 

.0199 

-.3283 

H44 

.3297 

.0408 

-.1325 

H45 

.2867 

.1249 

.0569 

H46 

.1248 

.1881 

.0505 

4 

Cll 

-.1412 

.2126 

.  2839 

C12 

-.2287 

.1915 

.  2007 

C13 

-.2002 

.1343 

.1036 

C14 

-.0952 

.1207 

.1267 

C15 

-.0587 

.1685 

.  2381 

5 

Hll 

-.1383 

.2519 

.3614 

H12 

-.2985 

.2132 

.  2091 

H13 

-.2463 

.1084 

.0312 

H14 

-.0538 

.0824 

.0736 

H15 

.0129 

.1710 

.2777 
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RESULTS 

The  molecular  geometry  and  the  numbering  system 
used  is  shown  in  Figure  16,  while  the  geometry  of  the 
manganese  atom  is  shown  in  Figure  17.  The  molecular 
packing  as  viewed  down  the  c  axis  is  depicted  in 
Figure  18.  The  drawings  were  made  using  the  programme 
ORTEP.  Table  XIX  gives  the  bond  lengths  and  angles 
within  the  molecule.  Some  intermolecular  contacts 
are  listed  in  Table  XX.  The  results  and  estimated 
errors  associated  with  them  were  calculated  using  the 
programme  0RFFE2. 
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Figure  16 


Molecular  structure  of  (tt-C^H  )  (OC) ~Mn (H) Si (CrH_ ) ~ 

o  £  o  5  J 
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Figure  17 


Geometry  around  manganese  atom 
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Figure  18 


Molecular  packing  projected  onto  ab  plane 
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Table  XIX 

o 

(a)  Bond  lengths  (A)  in 

(tt-CkH_)  (OC)  0Mn  (H)  Si  (CcH_)  - 
D  D  Z  O  D  J 


Mn-Si 

Si-H 

Mn-H 

Mn-Cl 

Mn-C2 

Cl-Ol 

C2-02 

Mn-Cll 

Mn-Cl2 

Mn-Cl 3 

Mn-Cl 4 

Mn-Cl 5 

ave  Mn-C 

C11-C12 

Si-C21 

Si-C31 

Si-C41 

Cl-H 


2.424  (2) 
1.76  (4) 
1.55(4) 
1.771(7) 
1.764(7) 
1.162(6) 
1.168  (6) 
2.140  (4) 
2.132  (4) 
2.136(4) 
2.146  (4) 
2.148  (4) 
2.140  (4) 
1.414  (3) 
1.883  (5) 
1.888  (5) 
1.886(5) 
2.08(4) 


C21-C22 

C21-C26 

C22-C23 

C23-C24 

C24-C25 

C25-C26 

C31-C32 

C31-C36 

C32-C33 

C33-C34 

C34-C35 

C35-C36 

C41-C42 

C41-C46 

C42-C43 

C43-C44 

C44-C45 

C45-C46 


1.393(7) 
1.398(7) 
1.385(8) 
1.365  (9) 
1.374  (9) 
1.384  (8) 

1.391  (7) 
1.396  (7) 
1.382(7) 
1.369  (8) 
1.382(8) 

1.392  (7) 
1.398(7) 
1.395  (7) 
1.415(9) 
1.378(9) 
1.353  (9) 
1.391(8) 


. 

' 

■ 
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Table  XIX 
(Contd . ) 


(b)  Bond  angles  (degrees)  in 
(tt-C5H5)  (OC)2Mn(H)Si(C6H5)3 


Mn-Hl-Si 

94(2) 

Mn-Si-Hl 

39(1) 

Si-Mn-Hl 

46(2) 

Cl-Mn-C2 

88.7(3) 

Si-Mn-Cp 

118.0(1) 

Cl-Mn-Cp 

122.6  (2) 

C2-Mn-Cp 

123.7(2) 

Hl-Mn-Cp 

122.0 (1) 

Cl-Mn-H 

77.0(1) 

C2-Mn-H 

109.0  (2) 

Mn-Si-C21 

108.1  (2) 

Mn-Si-C31 

114.1  (2) 

Mn-Si-C41 

115.9(2) 

C21-Si-C31 

106.6  (2) 

C21-Si-C41 

106.7(2) 

C31-Si-C41 

104.7  (2) 

01-Cl-Mn 

178.6  (6) 

02-C2-Mn 

176.3(5) 

C26-C21-C22 

116.3  (5) 

C21-C2  2-C23 

122.3(6) 

C22-C23-C24 

119.9(6) 

C23-C2  4-C25 

119.8(6) 

C24-C25-C26 

120.4  (6) 

C25-C26-C21 

121.4 (6) 

C36-C31-C32 

116.6(5) 

C31-C32-C33 

122.0  (5) 

C32-C33-C34 

120.4(5) 

C33-C34-C35 

119.5(5) 

C34-C35-C36 

119.9(5) 

C35-C36-C31 

121.6(5) 

C46-C41-C42 

116.2(5) 

C41-C42-C43 

121.8(6) 

C42-C43-C44 

118.1(6) 

C43-C44-C45 

122.3  (7) 

C44-C45-C46 

118.8  (6) 

C45-C46-C41 

122.9  (6) 

* 

'  * 

. 

, 
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Table  XX 

Intermolecular  contacts 


01-H33 

2.723  (5) 

01-H44 

2.829(5) 

01-H23 

3.288(5) 

02-H15 

2.937(4) 

02-H46 

3.270  (5) 

H22-H45 

2.646(1) 

H22-H11 

2.787(1) 

H25-H45 

2.662(1) 

H33-H44 

2.647(1) 

H34-H44 

2.610(2) 

H36-H15 

2.576(1) 

01-H24 

3.038(5) 

01-H35 

3.048(5) 

02-H25 

2.749  (4) 

02-H35 

3.049(4) 

02-H36 

3.284(4) 

H22-H46 

2.671(1) 

H23-H46 

2.849(1) 

H33-H34 

2.457(1) 

H34-H34 

2.415(1) 

H34-H13 

2.796 (1) 

H36-H11 

2.645(1) 
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DISCUSSION 

A  drawing  of  the  compound  (tt-C5H5)  (OC)2Mn(H) 
Si(C^H^)2  is  found  in  Figure  16.  The  geometry  about 
the  manganese  atom  (Figure  17)  can  be  described  as 
a  distorted  square  pyramid  and  that  of  the  silicon  atom 
as  a  distorted  tetrahedron.  The  molecule  has  a  hydrogen 

o 

cis  to  the  silicon  and  close  enough  (1.76 (4) A)  to 
interact  with  it. 

The  Si-C^.-H..)  distances  of  1.883(5)A,  1.888  (5)A 
6  5 
o 

and  1.886 (5) A  are  the  same  within  experimental  error. 

There  is  some  distortion  from  tetrahedral  geometry  about 

the  silicon  atom,  the  Mn-Si-C (C^H^)  angles  varying  from 

108.1(2)°  to  115.9(2)°  and  the  C (CrHc ) -Si-C (C,Hc ) 

6  5  6  5 

(different  rings)  from  104.7(2)  to  106.7(2).  Angles 
within  the  phenyl  rings  range  from  116.3(5)  at  the 
attached  carbons  to  122.3(6) .  It  will  be  noticed  that 
in  the  phenyl  rings  there  is  a  range  in  C-C  bond  lengths 

o 

from  1.365(9)  to  1.415 (9) A.  This  bond  length  range  and 
distribution  is  not  that  expected  for  a  single  popu¬ 
lation  of  formally  identical  bond  lengths  with  an 

o 

estimated  standard  deviation  of  0.008A  as  determined 
from  the  least  squares  refinement.  The  variation  of 
bond  lengths  follows  a  definite  pattern,  the  longest 
being  those  close  to  the  silicon-bonded  carbon  (average 

o 

1.395A).  The  shortest,  with  an  average  C-C  distance  of 


. 

' 
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o 

1 . 370A,  are  those  furthest  away.  A  similar  pattern 

5  3 

was  observed  by  Beauchamp  et  al .  and  results  from 
the  waggling  motion  of  the  rings.  The  apparent 
shortening  of  the  C-C  bond  is  due  to  the  same  effect 
as  discussed  in  the  hindered  rotor  section;  that  is, 
the  atoms  move  along  an  arc,  the  centre  of  which  lies 
near  the  silicon  atom.  Attempts  to  describe  the 
distribution  of  electron  density  resulting  from  such 
motion  by  means  of  anisotropic  temperature  factors 
would  give  an  error  in  the  centre  of  electron  density 
of  each  atom,  and  hence  in  the  calculated  bond  lengths. 
The  greatest  apparent  shortening  would  occur  between 
those  atoms  v/hich  have  the  greatest  motion,  i.e., 
CR3-CR4  and  CR4-CR5  where  R  is  the  number  of  the 
ring  (R  =  2,3,4). 

The  Mn-C  (tt-C^H^)  and  Mn-C(CO)  distances  were  not 

unusual  and  agree  with  those  of  such  compounds  as 

S4  24 

(C6H7)Mn(CO)  ^  4,  (C10Hg)Mn2  (CO)  g  ,  (tt-C^  )  Mn  (CO)  g 

and  [  (tt-C^H^CH^)  Mn  (CO)  2  ]  2CgH4  (As  (CH3)  2^  2  '  as  can 

be  seen  in  Table  XXI.  The  variations,  where  signi¬ 
ficant,  are  consistent  with  the  trends  that  would  be 
predicted  on  the  basis  of  the  relative  tt  acceptor 
properties  of  the  other  ligands  present.  In  all 
these  compounds  the  carbonyl  groups  lie  trans  to  a 
delocalized  tt  system.  One  would  expect  that  the 


* 


- 
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Mn-C(O)  distances  trans  to  other  carbonyls 

would  be  longer  than  those  trans  to  ir-ring  systems. 

This  more  dramatic  demonstration  of  the  trans  effect 

5  7 

can  be  seen  in  such  compounds  as  HMn(CO)c  , 

(CO)4Mn(H)  (P(C6H5)2)Mn<CO)  458,  Mn2(CO)lc|59,  and 

HRe^Mn (CO) in  which  Mn-C  distances  trans  to 

carbonyl  groups  are  longer  (Table  XXII) . 

The  Mn-Si  bond  length  of  2.424 (2) A  lies  between 

that  found  in  (OC)  4Mn  (Si  (C^)  2)  2Mn  (CO)  461  of  2.40A 

°  9 

and  that  of  2.497(5)A  in  Mn(Si(CHj0) (CO) _  .  It  is 

Jo  5 

significantly  shorter  than  the  calculated  single 

°  .  °  6  3 

bond  distance  of  2.56A  (using  1.17A  for  the  radius 

°  64 

of  silicon  and  1.39A  for  manganese),  and  would 
suggest  some  double  bond  character  between  the 
manganese  and  silicon. 

The  most  exciting  part  of  the  work  is  the 
location  of  the  bridging  hydrogen.  It  lies  in  an 
otherwise  empty  space  between  the  manganese  and  silicon 
atoms,  and  as  described  earlier  was  found  to  be  the 
highest  peak  on  all  low  angle  difference  maps.  It  is 

o 

bonded  to  the  manganese  with  a  distance  of  1.55  (4) A 

o 

or  ca  1.62A  when  corrected  for  the  use  of  a  spherical 

6  5 

electron  distribution  .  This  is  close  to  the  terminal 

o 

Mn-H  distance  of  1.60A  found  by  neutron  diffraction 
techniques  in  HMn(CO)  557,  1.60  (12)A  in  RhH  (CO)  (P  (C^)  3 )  34 


' 

- 
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Table  XXII 

Comparison  of  bond  lengths 


Molecule 

Mn-C 
trans  to 
carbonyl 

Mn-C 

not  trans  to 
carbonyl 

HMn (CO) _ 

D 

1.836(5) 

1.821(9) 

m2  (C0)  10 

1.830  (8) 

1.792(14) 

(CO)4Mn(H)  (P(C6H5)2Mn(CO)4 

1.840  (9) 

1.788  (9) 

HRe2Mn (CO) 14 

1.840  (15) 

1.80  (3) 

. 
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and  1.6(1)  in  CoH  (N0 )  (P  (C,-H  )  _ )  ^^  the  latter  two 

2  6  5  3  3 

being  done  using  X-ray  methods.  The  Si-H  distance 

o 

of  1.76 (4) A,  although  considerably  longer  than  the 

O  (Z  c, 

Si-H  bond  distance  of  1.48A  in  SiH^  ,  is  much 

shorter  than  would  be  expected  for  a  normal  Si-H 

intramolecular  contact.  Because  silicon  has  a  larger 

radius  than  carbon,  one  would  expect  the  C**»*H  contact 

to  be  less  than  that  of  Si****H.  This  is  not  the  case, 

Cl-H  being  2.08 (4) A  and  Si-H  being  1.76 (4) A.  This 

would  suggest  an  interaction  between  the  Si-H  and 

formation  of  a  weak  bond.  A  similar  type  of  inter- 

6  7 

action  was  found  in  bipy (OC) ^CZMoSnCH^CZ^  in  which 


the  Sn-CZ  distance  of  2. 805 (4) A  was  longer  than  would 
be  expected  for  a  bridging  chlorine  distance,  but  much 


shorter  than  Van  der  Waals  distance  (4.0A) 


The  bond 


bond 


length  evidence  would  suggest  a  bent  3  centre,  2  electron 
H  1 

^Si  wit*1  the  predominant  atomic  orbitals 
*  * 

lying  on  the  manganese  and  hydrogen  atoms. 

The  presence  of  such  a  bridge  would  help  to  explain 
the  different  spectral  and  kinetic  properties  observed 


between  this  compound  and  the  trihalogenated  silyl 

derivatives.  In  the  structure  of  ( tt-C^^)  (OC)  Fe  (H) 

6  8 

(SiCZ3)2  ,  the  hydrogen  atom  was  not  located.  However, 

o 

if  an  Fe-H  distance  of  ca  1.6A  is  assumed,  the  hydrogen 


atom  lying  midway  between  the  Si's  and  ca  120  from 


. 


. 
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Fe-Centre  (ring) ,  then  one  would  obtain  a  Si-H 

o 

distance  of  2.3A.  This  is  greater  than  the  Cl****H 

o 

distance  of  2.08 (4) A  and  implies  that  there  is  no 

Fe^  ^Si  brid9e*  The  more  electronegative 
chlorines  pull  sufficient  electron  density  from  the 
silicon  to  destroy  any  Si-H  interaction.  Alter¬ 
natively,  taking  the  electrostatic  view  of  the  hydrogen 
6  —  64-  6  — 

bond  X  -  H* • • *Y,  the  requirement  is  that  the  Y  group 
(in  this  case  silicon)  have  a  partial  negative  charge. 
In  SiCl^,  the  highly  electronegative  chlorine  would  be 
expected  to  produce  a  slight  positive  charge  on  the 
silicon,  which  while  improving  the  electron  acceptor 
properties  as  far  as  transition  metal  interactions 
were  concerned,  would  reduce  the  ability  to  hydrogen 
bond.  It  would  thus  appear  that  the  difference  in 
physical  properties  of  the  two  types  of  compound  is 
due,  in  part  at  least,  to  the  presence  or  lack  of  a 
bridging  hydrogen.  The  chlorinated  silyl  compounds 
have  a  terminal  M-H  bond  and  clearly  show  this  in  the 
infrared.  The  phenyl  substituted  silyl  class  have  a 
Mn-H  bond  that  is  close  to,  but  not  exactly  the  same 
as,  a  terminal  bond.  This  difference,  due  to  the 
Si • • • *H  interaction,  could  help  to  explain  the  lack 
of  a  peak  in  the  infrared  M-H  terminal  region,  but  the 
appearance  of  a  weaker  peak  in  the  Raman. 


■ 

' 
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69 

The  preliminary  results  of  a  kinetic  study  of 
the  displacement  of  a  silane  by  a  phosphine,  which  may 
be  represented  by  the  following  general  chemical 
equation: 

nCcH_Mn  (CO)  „SiHR0+P  (C,HJ  _  ->  TTCcHirMn  (CO )  0P  (C  CH_  )  _  +  SiHR_ 

O  D  Z  D  ODD  DO  Z  O  D  D  D 

may  be  rationalized  in  terras  of  the  structural  results 
reported  here.  In  the  case  of  the  triphenyl  silyl 
derivative,  the  dominant  reaction  proceeds  via  an  initial 
dissociation  in  which  triphenyl  silane  is  eliminated  and 
the  second  step  involves  very  rapid  phosphine  attack  on 
the  electron  deficient  species  ttC^H  Mn  (CO)  ^ .  The 
reaction  is  the  first  order  in  transition  metal  complex 
and  zero  order  with  respect  to  phosphine.  However,  for 
the  trichlorosily 1  derivative  the  main  reaction  is  found 
to  be  bimolecular.  The  formation  of  a  weak  silicon 
hydrogen  bond  corresponds  to  a  partial  activation  of 
the  leaving  group  and  significant  formation  of  the 
silicon  hydrogen  bond  in  the  trichlorosilyl  derivative 
occurs  only  with  the  involvement  of  the  entering  phos¬ 
phine  ligand. 


. 


- 


Chapter  IV 


THE  CRYSTAL  STRUCTURE  OF  DIHYDROOCTACARBONYL  DIRHENIUM 
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INTRODUCTION 

There  has  been  great  interest  lately  in  the  bonding 

and  geometry  of  M-H-M  bridges.  The  presence  of  the 

hydrogen  should  result  in  a  lengthening  of  the  M-M 

distance,  since  a  2  electron/3  centre  bond  is  postulated. 

6  0 

This  was  observed  in  the  structure  of  HRe2Mn (CO) 

o 

in  which  a  Re-Re  distance  of  3.39A  was  found.  This 

o  o 

is  longer  than  the  Re-Re  distance  of  2.96A  to  3.02A 

found  in  such  structures  as  Re^ (CO) ^  .  A  similar 

result  was  found  by  Dahl  in  the  structure  of 

—  71 

(HCr2  (CO) ^q)  .  In  this  molecule  a  Cr-Cr  distance 

o 

of  3.41A  was  obtained.  If  a  linear,  symmetrical 
M-H-M  bond  is  assumed,  both  compounds  have  a  M-H 

o 

distance  of  1.70A  which  is  similar  to  the  neutron 

°  .  =72 

diffraction  result  of  1.68 (1)A  for  Re-H  m  ReH^ 

It  should  also  be  possible  to  have  a  bent  M-H-M  bridge. 
One  would  then  expect  a  M-M  bond  shorter  than  that 
found  in  the  linear  bridge,  but  longer  than  that  in 

the  unbridged  M-M  bond.  This  has  been  found  in  such 

-  73 

structures  as  (H2Re3 (CO) 12  )  in  which  bridged  Re-Re 

o  o 

distances  of  3.181 (7)A  and  3.173 (7)A  were  found,  and 

in  (H,-Re .  (CO)  ,  ~)  7  having  an  average  Re-Re  distance 
6  4  12 

o 

of  3.160 (7) A.  Having  established  a  pattern  for 
single  bridged  metal-metal  bonds,  one  wonders  what 


£ ■  « 3  1  C-  .  i  ,ic  th 
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would  occur  in  doubly  hydrogen  bridged  bonds  such 

75 

as  occur  in  H2Re2(CO)g  and  H2OSg (CO)  .  These 
bridges  are  now  4  centre  rather  than  3  centre 
systems.  Bonding  schemes  such  as  protonated  metal- 
metal  double  bonds  (M=M)  have  been  suggested  as  well 
as  a  4  centre  delocalization  of  electrons.  In  an 
attempt  to  understand  the  bonding  in  such  systems, 
the  crystal  structure  of  H2Re2 (CO) g  was  undertaken. 


Ill 


EXPERIMENTAL 

The  compound  HoRe(C0)o,  molecular  weight  598.4 

2.  o 

gm/mole,  was  prepared  by  J.K.  Hoyano  and  was  re¬ 
crystallized  from  a  mixture  of  dichloromethane  and 
hexane  to  obtain  crystals  suitable  for  single 
crystal  X-ray  diffraction  studies.  The  yellow 
crystals  were  mounted  parallel  to  the  elongated  edge 
of  the  crystal.  Rotation,  Weissenberg  (0k£,  h0£- 
h7£,  hkO-hkl)  and  precession  (h0£,  0k£)  photographs 
showed  the  crystal  to  be  monoclinic  with  systematic 
absences  h(H :  l  =  2n+l  and  OkO :  k  =  2n+l,  fixing 
the  space  group  at  P2^/c.  The  unit  cell  was  found 
to  be  a  =  8.96 (2 ) A,  b  =  11.62(2)A,  and  c  =  12.85(2)A, 

3  =  109.2(1)°  from  precession  and  Weissenberg  photo- 
7 

graphs  .  The  density  was  determined  experimentally 
by  means  of  flotation  in  a  solution  of  Clerici's 
Reagent  and  water  to  be  3.16  gm/cc  which  agreed  with 
the  calculated  density  of  3.18(2)  gm/cc  for  four  mole¬ 
cules  per  unit  cell.  This  would  require  that  the  mole¬ 
cules  lie  in  general  positions  in  the  unit  cell.  Data 

o 

within  the  1.0A  sphere  were  collected  on  the  PAILRED 
automatic  diffractometer  using  Weissenberg  geometry 
and  a  moving  crystal/stationary  counter  technique 
(a)  scan)  .  The  hexagonal  crystal  was  mounted  along 
the  b  axis,  and  13  layers  (h0£-hl2£)  were  collected 


. 
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using  MoK  -radiation .  The  X-ray  beam  was  monochromated 

using  a  graphite  single  crystal  (002  reflection) . 

A  scintillation  counter,  aperture  2.5°,  was  used,  with 

the  pulse  height  analyser  centred  on  the  MoK  peak  and 

set  to  accept  90%  of  the  peak.  A  half  scan  of  1.5° 

for  the  inner  reflections  (sin  0/A  ,<0.20)  and  1.0° 

for  the  outer  reflections  was  used  at  a  scan  speed  of 

2.5°/min.  The  scan  widths  were  increased  to  2.0° 

and  1.5°  for  inner  and  outer  reflections  for  the  higher 

levels.  The  peak  was  scanned  three  times  or  until 

the  peak  count  passed  4,000  counts.  A  stationary 

background  was  counted  for  twenty  seconds  on  each 

side  of  the  scan.  Because  the  crystal  was  found  to 

decompose  in  light,  and  to  a  lesser  extent  in  radiation, 

the  data  were  collected  in  a  darkened  room.  Twelve 

standard  reflections  were  measured  at  the  end  of  each 

layer  and  a  decomposition  correction  applied  as  a 

function  of  time  and  sin  0/A.  The  intensities  were 

calculated  using  I  =  P-t(B^+B2)  where  P  is  the  peak 

count,  B^  B2  the  background  count,  and  t  is  the  ratio 

of  peak  count  time  to  that  of  total  background.  Four 

hundred  and  twenty-two  reflections  were  rejected  using 

the  criteria  (1)  1=0  and  (2)  I  <  3AI  where  AI  = 

(P+t2  (B  +B  ))i$.  Absorption  corrections  were 
JL  2 

applied  (y (Mo)  =  217. cm-1,  max  crystal  dimensions: 

.14  mm  x  .11  mm  x  .14  mm,  mean  yr  =  1.4). 


•  .1  • 
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SOLUTION  AND  REFINEMENT 


The  positions  of  the  rhenium  atoms  were  obtained 

from  a  three  dimensional  Patterson  map  as  Rel: 

x  =  0.295,  y  =  0.130,  z  =  0.273  and  Re2:  x  =  -0.046, 

y  =  0.130,  z  =  0.180.  A  structure  factor  calculation 

using  the  form  factors  calculated  from  Cromer's 
.  .  50 

coefficients  gave  reliability  factors  of  R^  = 

27.9%,  R2  =  31.3%,  where  R^  =  z||F0|-|Fc||  • 

E  I  F  I 
1  o  1 

r  \ 


r2  - 


IW(|F0|-|FC|) 


Ew  F 


o 


^  ,  w  =  l/a(F)2  ancj  a  (F) 


(P+t^CB^+B^)  +  P^I)  2,  p  =  0.06.  A  difference  map 
phased  by  the  two  rhenium  atoms  gave  the  positions  of 
all  other  non-hydrogen  atoms.  Four  cycles  of  full 
matrix  least  squares  refinement,  minimizing  the 

ii  i  2 

function  Ew ( | F  | — | F  |)  with  all  atoms  isotropic, 

o  c 

reduced  the  R  factors  to  R^  =  16.0%  and  R2  =  20.4%. 

It  was  noticed  that  four  layers  (hl£,  h2£,  h3£  and 
h9£)  had  much  higher  R  factors  (of  the  order  of  20%) 
than  the  others  (-10%).  For  this  reason,  data  for 
these  layers  were  recollected  using  a  second  crystal 
of  similar  dimensions  (.12  mm  x  .13  mm  x  .10  mm). 

The  data  were  treated  as  for  the  first  crystal, 
corrected  for  absorption  and  combined  with  the  previous 


* 
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data  using  a  second  scale  factor.  Inclusion  of  the 
anomalous  dispersion  corrections  for  rhenium  (Af'  = 
“1.40,  Af"  =  7.70)  in  the  calculated  structure  factors 
and  further  refinement  dropped  the  reliability  factors 
to  =  10.3%,  R2  =  12.0%.  A  difference  map  at  this 
point  showed  peaks  of  3-4  electrons  near  each  of  the 
rheniums,  suggesting  an  error  in  the  absorption 
correction.  Since  any  anisotropic  refinement  would 
be  meaningless  for  this  reason,  and  also  because  of  the 
large  decomposition  correction  (-30%)  applied,  refine¬ 
ment  was  terminated  at  this  point.  In  the  final 
cycle,  all  atoms  isotropic,  no  shift  was  greater  than 
1/10  the  corresponding  e.s.d.  The  observed  and 
calculated  structure  factor  amplitudes  (xlO)  are 
listed  in  Table  XXIII.  The  positional  parameters 
and  temperature  factors  are  given  in  Table  XXIV. 
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Table  XXIII 

Observed  and  calculated  structure  factor 
amplitudes  (xlO)  in  units  of  electrons 
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4391 

l 

1 

2496 

2863 

0 

0 

3766 

3840 

10 

-8 

538 

597 

9 

-6 

817 

882 

0 

2 

918 

802 

7 

4 

491 

837 

1C 

-7 

550 

386 

0 

8 

2753 

3144 

1 

2 

1  794 

2087 

0 

2 

1787 

1650 

10 

-9 

1318 

1311 

9 

-7 

ll  35 

1235 

0 

4 

397 

344 

7 

5 

1038 

940 

10 

-8 

806 

768 

0 

10 

1487 

16  04 

1 

3 

2033 

2486 

0 

8 

907 

916 

to 

-10 

1344 

14  75 

9 

-ll 

5  52 

474 

0 

6 

473 

541 

7 

15 

0 

612 

ll 

15 

0 

1624 

0 

-2 

1843 

1948 

l 

4 

2001 

2302 

0 

1C 

1516 

1606 

10 

-ll 

1239 

1354 

9 

-12 

601 

457 

0 

8 

1316 

1328 

7 

-1 

164? 

1523 

1 1 

-1 

1123 

965 

0 

-4 

588 

572 

l 

3 

1419 

1617 

0 

1? 

1360 

1404 

10 

-12 

1144 

1193 

10 

3 

640 

655 

0 

-2 

l  COT 

437 

7 

-3 

2164 

196  1 

1 1 

-3 

473 

312 

0 

-6 

576 

657 

1 

6 

92  7 

1052 

c 

-2 

3344 

3031 

10 

-13 

736 

817 

10 

5 

963 

94  7 

0 

-4 

1  123 

944 

7 

-5 

20)4 

1744 

11 

-9 

1112 

1008 

0 

-8 

786 

721 

1 

7 

601 

630 

c 

-4 

2482 

2455 

11 

0 

1255 

11  16 

10 

15 

0 

4  75 

n 

-4 

1491 

1  395 

7 

-  7 

1052 

1011 

0 

-10 

446 

477 

l 

-1 

1957 

1557 

0 

-6 

710 
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11 

1 

1030 

935 

10 

-l 

1275 

1055 

0 
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2566 
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1 
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0 
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1143 

1228 

l 
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1234 

1145 

0 

-8 

1764 

1  783 

1 1 

2 

1135 

1150 

10 

-2 

867 

753 

0 

-12 

182  7 

1903 

n 

15 

r. 

76  4 

0 

-14 

1639 

1819 

l 

-4 

1470 

1449 

0 

-10 

1944 

1991 

11 

3 

524 

625 

10 

-3 

1903 

1728 

0 

-  1 4 

1147 

ll®6 

8 

-l 

493 

792 

n 

54  7 

391 

1 

2 

444 

526 

l 

-5 

l  75C 

1860 

0 

-12 

950 

1089 

11 

15 

0 

2965 

10 

-4 

1061 

1739 

0 

2904 

7440 

a 

-3 

582 

398 

0 

-4 

956 

1025 

1 

7 

789 

912 

l 

-6 

1346 

1437 

l 

0 

289  1 

2686 

11 

-1 

4402 

3622 
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-5 

1466 
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24  14 

2148 
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582 
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1 

B 
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1216 
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-7 

1310 

1413 

1 

1 

2550 

2501 

11 

-3 

1925 

1907 

10 

-6 
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2 

l  785 

1556 

9 

C 
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-0 

1069 

993 

1 

4 

1007 

1288 

l 
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565 
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1 

2 

2068 

2097 

11 

-7 

1505 

1557 

1C 

-7 

756 
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) 

901 

940 

9 

-l 

704 

716 

.  j 

656 
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1 

10 

46  3 

1043 

2 

1 

731 
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1 
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1643 

1603 

11 

-9 

1589 

1708 

10 
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701 
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- 1 

31  <.5 

2552 

9 

-? 

2505 

2105 

.j 

881 
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1 

11 

65  3 
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2 

2 

1721 

1609 

1 

4 

950 

887 

11 

-ll 

1032 

9  70 

1C 

-ll 

813 

831 

-2 

2330 

1975 

9 

-4 

3334 
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.) 

1120 
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1 

-3 

2489 
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3 

1904 

2162 
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1032 

992 
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950 
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2392 

2302 

2 

7 

1601 

1706 

1 

9 

11(3 

1211 

11 
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Table  XXIV 

Positional  parameters  for  H  Re  (C0)o 

Z  2  8 


Atom 

X 

y 

z 

Isotropic 

B 

Rel 

.2948(2) 

.1289(2) 

.2691(1) 

1.69  (5) 

Re2 

-.0467(2) 

.1190 (2) 

.  1800  (2) 

1.80(5) 

Cl 

.  2948  (78) 

. 2262  (53) 

.1432(50) 

2.6(11) 

01 

.3114 (65) 

.2936(42) 

.0798(41) 

4.0(10) 

C2 

.4701(67) 

.0340(48) 

.2412 (43) 

1.8(9) 

02 

.5621(65) 

-.  0183  (40) 

.  2354  (40) 

3.8(9) 

C3 

.4502 (65) 

.  2249  (42) 

.3710(43) 

1.7(9) 

03 

.5460  (64) 

.  2812  (43) 

.4340(41) 

4.0(10) 

C4 

.  2994  (71) 

.0220  (47) 

.  3921  (46) 

2.1(10) 

04 

.3135(69) 

-.0325(48) 

.4675(46) 

2.0(10) 

C5 

-.0512(85) 

.  0178  (56) 

.3111(55) 

3.0(13) 

05 

-.0721(61) 

-.0365  (40) 

.3680 (41) 

3.7(9) 

C6 

-.2190(88) 

.2035(56) 

.2005(54) 

3.0(12) 

06 

-.3072(59) 

.2611(40) 

.2199 (38) 

3.4(9) 

C7 

-.1934 (95) 

.0163 (61) 

.  0827  (58) 

3.6(13) 

07 

-.2868(76) 

-.0499 (51) 

.  0298  (48) 

5.3  (12) 

C8 

-.0451(73) 

.2144  (51) 

.0457(48) 

2.3 (10) 

08 

-.0743  (73) 

.  2576  (50) 

-.0362(49) 

5.1(12) 

• 
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RESULTS 

A  drawing  of  the  molecule  is  shown  in  Figure  19, 
giving  the  numbering  scheme  used.  The  diagram  was 
drawn  using  the  Calcomp  plotter  and  Carroll  K.  Johnson's 
programme  ORTEP.  The  bond  lengths  and  angles  within 
the  molecule  are  found  in  Table  XXV,  while  Table 
XXVI  gives  some  intermolecular  contacts.  These 
distances  and  their  associated  errors  were  calculated 
using  0RFFE2. 
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Figure  19 


View  of  H2Re2(CO)g  down 


axis 
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Table  XXV 


(a)  Bond  lengths  (A)  in  H0Re„ (CO) 0 

A  A  o 


Re  1- Re 2 

2.896(3) 

Cl-01 

1.17  (7) 

Rel-Cl 

1.97  (6) 

C2-02 

1.05(7) 

Rel-C2 

2.05(6) 

C3-03 

1.17(7) 

Rel-C3 

1.92(5) 

C4-04 

1.13(8) 

Rel-C4 

2.00(6) 

C5-05 

1.03(7) 

Re2-C5 

2.06(7) 

C6-06 

1.13(8) 

Re2-C6 

1.92(7) 

C7-07 

1.17  (9) 

Re2-C7 

1.91(7) 

C8-08 

1.12  (7) 

Re2-C8 

2.06(6) 

(b)  Bond  angles  in 

H2Re2 (CO) 

8 

-Rel-C2 

87.  (2) 

C5-Re2-C6 

87.  (3) 

-Rel-C4 

176.  (2) 

C5-Re2-C8 

178.  (2) 

-Rel-C3 

90.  (2) 

C6-Re2-C7 

90.  (3) 

-Rel-C4 

90.  (2) 

C7-Re2-C8 

89.  (3) 

-Rel-Re2 

89.  (2) 

C5-Re2-Rel  90.(2) 

(c)  Intramolecular 

non-bond 

contacts 

C1-C8 

2.89(9) 

C4-C5 

2.97  (10) 

01-08 

3.31(8) 

04-05 

3.27(8) 

C1-C2 

2.78(8) 

C5-C6 

2.74  (9) 

C2-C3 

2.82(8) 

C6-C7 

2.70  (10) 

C3-C4 

2.77  (7) 

C7-C8 

2.78(9) 

' 
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Table  XXVI 

o 

Intermolecular  contacts  (A) 


Cl-07 

3.40(8) 

Cl-02 

3.40  (8) 

C1-C7 

3.93(9) 

01-04 

3.13 (8) 

01-07 

3.14(8) 

01-05 

3.15(7) 

01-02 

3.15(7) 

01-08 

3.31(8) 

01-C8 

3.21(8) 

01-C3 

3.32(8) 

01-06 

3.32(7) 

C3-04 

3.30(7) 

C3-06 

3.38  (7) 

C3-07 

3.44(8) 

C4-05 

3.31(8) 

C4-04 

3.34(8) 

C4-06 

3.37  (7) 

C5-08 

3.31(9) 

C6-08 

3.23  (9) 

C2-07 

3.33  (8) 

C2-06 

3.37(7) 

C2-C6 

3.59  (9) 

C2-04 

3.60  (8) 

02-03 

3.12  (7) 

02-05 

3.16(8) 

02-C3 

3.27(7) 

02-C5 

3.30(9) 

02-C6 

3.36(8) 

02-07 

3.36(8) 

02-C7 

3.41(9) 

03-07 

3.19  (8) 

03-04 

3.24(8) 

03-08 

3.33(8) 

03-06 

3.43  (7) 

04-04 

3.26(12) 

04-06 

3.38(7) 

05-C8 

3.15  (7) 

06-08 

3.15(7) 

. 


■ 

. 
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DISCUSSION 

The  compound  Re^H^ (CO) g  shows  octahedral  geometry 
about  the  rhenium  atoms  if  one  assumes  that  the  hydrogen 
atoms  occupy  the  empty  region  between  the  Re  atoms  and 
complete  the  octahedron.  The  molecule  as  a  whole  has 
symmetry,  although  none  is  required  crystallo- 
graphically . 

There  is  little  deviation  from  90°  of  the  C-Re-C 

angles,  the  largest  being  -3.(2)  degrees.  The  Re-C 

°  ° 

bonds  range  in  length  from  1.91A  to  2.06A.  Those 
being  cis  to  the  hydrogens  are  longer  than  those  trans 

o 

(average  2.02 (3) A  as  compared  to  1.95(4)).  Since  the 
hydrogen-rhenium  bond  is  weaker  than  rhenium-carbon, 
the  bonds  trans  to  the  hydrogen  would  be  stronger  than 
those  trans  to  carbonyl  and  thus  have  shorter  bond 
lengths.  The  hydrogen  atoms  are  assumed  to  occupy 
the  two  remaining  octrahedral  sites.  If  there  is 
undistorted  octahedral  symmetry  (i.e.,  H-Re-H  of  90°), 

o 

then  a  Re-Re  distance  of  2.896 (3) A  would  require  a 

o 

Re-H  distance  of  2.05A.  This  is  rather  longer  than 

=  72 

the  average  Re-H  distance  (terminal)  in  ReHg  of 

o 

1.68(1) A  from  neutron  diffraction  work,  and  would 
suggest  that  there  is  a  slight  distortion  from  octa¬ 
hedral  symmetry  by  the  reduction  of  the  H-Re-H  angle 


. 

. 

- 
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in  order  to  shorten  the  Re-H  distance.  The  Re-H 

o 

bond  length  lies  between  two  extremes  -  2.05A  for 

o 

undistorted  octahedral  symmetry,  and  1.45A  if  the 
M-H-M  unit  is  linear.  Most  likely  it  is  close  to 

o 

those  found  in  other  Re-H  structures  of  -1.70A.  This 
would  require  a  H-Re-H  angle  of  63°.  The  Re-Re 

o 

distance  of  2.896 (3) A  is  shorter  than  would  be 
expected  for  a  Re-Re  unbridged  bond.  A  Re-Re  bond 

of  3.035A  was  found  in  (^2^63  (CO)  2.96,  2.96 

0  =70  0 

and  3.02A  in  Re.(C0)1c  and  3.000A  in 

4  16 

(Si (C,Hc) 0) 0Re0 (CO) q77 .  A  similar  shortening  of 
6  d  2  z  z  o 

75 

the  metal-metal  bond  was  observed  in  H20s2 (CO) g  , 

also  containing  a  doubly  hydrogen  bridged  Os-Os  bond. 

This  shortening  of  M-M  bonds  is  the  opposite  to  the 

effect  of  bond  lengthening  found  in  singly  bridged 

bonds,  and  suggests  that  the  bonding  of  the  M\H/M 

H 

unit  is  very  different  from  that  of  the  \M 

group . 

In  order  to  explain  the  bonding  of  the  rhenium 

and  hydrogen  atoms,  a  molecular  orbital  calculation 

was  done.  If  one  considers  only  the  rhenium-hydrogen  unit, 

7  7 

then  the  following  orbitals  can  be  drawn  ,  assuming 
no  interaction  with  the  carbonyl  groups: 


' 


■ 

' 
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Re  1 


->■  x 


It  should  be  borne  in  mind  that  such  an  assumption, 

as  well  as  the  choice  of  orbitals,  can  prejudice  the 

results  and  may  invalidate  them.  Using  symmetry, 

42 

the  following  orbitals  are  permitted  : 


A 

g 


+  B 


3u 


'  =  A  +  B0 

H  g  2u 


The  orbitals 


*1 

*1 

*2 

*3 

*4 

*5 

*6 


rxy 

=  Bn  +  B0 
lg  2u 

to  <|>  g  are 

defined  as  : 

follows 

5d  0 

0  orbital 

on 

Rel 

X 

y 

5d  0 

0  orbital 

z 

on 

Re  2 

X 

y 

Is 

orbital 

on 

HI 

Is 

orbital 

on 

H2 

5d 

orbital 

on 

Rel 

xy 

5d 

orbital 

on 

Re  2 

xy 

t 


■ 
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Then  the  A  orbitals  can 

g 

and  ^  " 

Symmetry  considerations  ; 
•  • 

Using  Huckel  molecular  o 
equations  may  be  set  up 

ha'a,_e  ha'a" 
ha'a"  ha"a"~e 

Defining  a 

» 

V  =  't’5fi°*5  “ 

< 

1 

» 

=  V°*6  6 

a 

Br.  =  V°*6  6R- 

« 

SR'H  ‘  1  V  +  S 

4 5^° ^ 3  ^H- 


be  expressed  as 


it  that  these  may  combine. 
7  8 

al  theory  ,  the  Slater 
the  determinant  obtained: 


<!>]_ H ° 

-L  3  t  = 

-e- 

N) 

o 

+  28t 

*3H° 

*33t  = 

uH° 

<(■  43  T 

4>r » ° 

= 

^  J 

*2r 

0-j_  3  T 

4^° 

h3T  =  j 

<j)43T 

o 

CCS 

CN 

-e- 

4>_8t  = 

0  j 

* 

hh° 

^4  3  T 

-e- 

U) 

tcj 

o 

4>43t  = 

• 

ua° 

■t-33T 
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the  determinant  can  then  be  rewritten  as 


aR  +  3Rr“E  23 


R-H 


23R-H  aH  +  3HH  E 


=  0 


If  the  assumption  is  made  that  the  overlap  between  the 
rhenium  and  hydrogen  atoms  is  very  small,  then  3R_H  “  0 
and  the  expression  simplifies  to 


»R  +  eRR"E  ° 


0 


“H  +  PHirE 


0 


which  can  be  very  easily  solved  to  give 

Ei  = « 
and  E^  ~  a 


R 

H 


+  3 


RR 


+ 


Similar  calculations  can  be  carried  out  for  the  other 
set  of  orbitals,  s 


ijj '  B 


2u 


=  1 

/2 


*5  +  't’e 


and  V'B2U  =  1_ 

/2 


\ 


The  energies  were  calculated  to  be 


-b  + 


/b2-4c 

2c 


and  E,  =  -b  -  Ad2-4c 
q  2c 


X 
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where  b  =  aRl  +  3R,  +  aR  -  3R 

and  c  =  46 2R, H  +  aRlaH  -  aR,6H  +  aH^R> 

2 

If  the  assumption  is  made  that  the  3  and  cross 
are  much  smaller  than  other  terms,  then  and 
approximated  to 

E3  “  °H  ~  SH  +  BR'H 
and  E4  =  aRl  +  6Rl  -  6R,H 

There  is  only  one  orbital  combination  with 


* 


B 


3u 


=  1 

/2 


*1  ^2 


E5  aR  3RR 


and  one  with  B^  symmetry: 


B 


=  1 
lg  /2 


^5  “  *6 


E6  aR  0R'R' 


3h3R' 

3  terms 

E.  can  be 

4 


symmetry : 


- 
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The  following  energy  level  diagram  will  be  obtained: 


B 


3u 


aR  eR  B3u 


aR  3R 


B 


Ig 


aR'  6R' 


B 


B 


Ig 


2u 


aR'  3R' 


aR'  +  3R'  3R'H 


a 


R' 


B 


2u 


aR-+  SR' 


B 


2u 


6 


H 


B 


2u 


3 


R 


+  3 


R'H 


aH  + 


H 


aH  + 


3 


H 


where  the  order  of  the  levels  B^  and  the  higher  B2u  is 
unknown  and  depends  on  the  relative  sizes  °f  V  and 
3R|R.  If  3r,h  is  greater  than  2BrI ,  then  Blg  will  lie 
below  B2u.  There  are  eight  electrons  which  occupy 
these  orbitals  so  that  the  bottom  four  will  be  filled. 

The  bonding  of  the  Re-H  unit  will  be  (1)  both  Ag  orbitals 
filled  (that  is,  (a)  the  electrons  localized  on  the 


•  I 

. 
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hydrogen  atoms,  and  (b)  a  Re-Re  a  bond),  (2)  the  bonding 

orbital  filled  in  which  the  electrons  are  localized 
2u 

over  both  rheniums  and  both  hydrogens,  forming  Re-H 
bonds : 


and  (3)  one  of  the  antibonding  orbitals  depending  on 


the  relative  size  of  3  ,  and  3  ,  . 

K  K  n 


If  Bn  lies  below 

ig 


B2u'  the  electrons  occupy  a  Re-Re  antibonding 

orbital : 


\  / 

Rel  Re2 


/  w  \ 


Otherwise,  the  two  remaining  electrons  will  occupy  a 
Re-H  antibonding  orbital  in  which  the  electrons  are 
delocalized  over  all  four  atoms: 
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The  bonding  can  thus  be  considered  as  a  Re-Re  single 

bond  with  delocalization  of  the  remaining  electrons 

over  both  the  rheniums  and  hydrogens,  rather  than  a 

protonation  of  a  Re=Re  double  bond.  Examination  of 

the  energy  level  diagram  shows  that  the  hydrogens  act 

as  a  stabilizing  unit  even  if  there  were  no  Re-H  bond. 

++ 

Removal  of  the  hydrogens  to  give  Re-Re  would  result 
in  no  bond  at  all,  since  all  four  of  the  remaining 
levels  would  be  filled. 
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SUBROUTINE  ROTOR 

COMMON  TITLE!  18)  ,MDDE,  I ‘IV, NAD,  I  SA  N  ,  NP'J  ,  NEW  (  5)  , NESQ,  NFG'JR, 
INCUR , KARD, I D, LV , NLV , MS , MU , MU,  ML  SC ,KR I T , RP  AR , NATE,  I CENT 
2  ♦  N RIG,  NON R  ,  NH  I  M 

common  a  ,  b  ,  c ,  al  ,  b e  ,  ga  ,  a q  ,  bq  ,  c o  ,  alq,  3EQ,GAQ,NSL,SUMDL 

COMMON  NP  AR  ,  N V AR  ,  NC  NT  , KS  E  L ( 02 1  ) 

CO  MM  ON  S(  50  )  ,  BO  ,  GP  (  70  )  ,  GQ  (  60  )  ,  GP.  (  60  )  ,MF  ( 60) , MG (60 ) , MI (60) , 

2X ( 140 ) , Y{ 140) ,7 ( 140) , B  I  1  (  LAO) , B2? ( 60) , R33 ( 60)  ,B1  2(  60)  , 

2B 13 { 60)  ,  P23 (60)  , NA ( 70 )  , NAMF ( 140 )  , X C ( 1 0 )  ,  YC(  10 ) , ZC { 10)  , 

3P(  10)  ,R  (10)  ,U( 10 ) , BPB ( 10, 3) ,NS  TEP ( 10) , NURD (  10)  ,  R I D (  10) 

C  r  M M  UN  ML  S  W  ,  AW  ,  BW  ,  C  W ,  DW  ,  E  W ,  F W ,  G U  ,  HW  ,  S  I  Z  E  1  ,  S  I  Z  E  2 
COMMON  SAEP.  (10,10)  ,  NS  AFP  (  10, ID  )  ,  S  A  E  S  (  10)  ,  N  S  A  F  S  (  1 0  )  ,  S  A  E  F  (  10)  , 
lNSAFF(lO) 

COMMON  NTRANK,  12,12)  ,NVTN(4,36) 

COMMON  NF,FM(  32,20)  ,F(20)  ,FP{ 2  0)  ,FPP(20)  , RHO, RHOSQ 

COMMON  FOBS ,  »TAl , S FC , AOBS , AC A _ , BOBS , BC AL , OF L , S I GMA , EXT i , E X T2 

1  t MSI , MR  E J , Q , NH , MK , ML, TH , TK, TL ,DA(18Q), 08(130), DF (180) 

COMMON  C.HG(  921)  ,  STD  (921  )  ,  EPS(  1  80)  ,0  I  AG  (  130)  ,VU80) 

COMMON  All, A12, A  13, A2  2, A 23, A3  3 , BA  1 , 8A2 , 8 A  3 , B A 4 , B A5 , BA6 
COMMON  SN ( 50 ) ,BGN,GPN(70 ) ,GQN( 60) , 3RN(60) ,XN( 140) ,YN( 140) , 

1ZM140  )  ,  B  11N  (  140  )  ,  B  2  2 N  (  60)  ,  B33M  (  60  )  ,  3 1 2N  (  60 )  ,  B  1 3N  (  60  )  , 

2  B  2  3  N ( 6  0  )  ,XCN(  10) ,YCM( 10), ZCN(IO) , PN( 10) ,RN(  10)  ,UN( 10) , 

3BP  B ( 10, 8) , R I D N ( 10) 

COMMON/ SPSY/I, J,N,RH,RK,PL,T,RHH,RKK,RLL,RHK,RHL,RKL,NR 
COMMON  /PDC/SPK,SPL,SPP , SPR , SPU , CPH  ,  CPK , CPL , CPP  ,CPR , CPU , WL » 

1  SACAL,  WB, SBCAL,K,L,TPI , JF,LC  ,SPH, INV1 
C  OMMGN/BFSS /OBJ  B ( 1O,20),LM(6),BR,BEA,TE,TD,BDDtBD,BIA(1O,2O) 
1,DBI3 ( 10) ,BIO( 10) , JPA 
DIMENSION  WL(20) ,WB(20) 

FQUIVAL ENCE (B I , Bll ) , ( BIN, B11N) , ( AR, SN) 

DIMENSION  AR (20100) 

0  r  MENS  I  ON  8  T ( 1 40 ) , BI M ( 140 ) 

DIMENSION  SUM ( 6 ) , DELMS ( 2,6) , OE  L  MC ( 2 , 6) , DF LC { 2 , 3 ) , OE  L  S ( 2,3) 

REAL  MTERM 

RFAL*8  CC,DD,AP,AP2 

COMMON/ RIGID/CLIMIT 

TE  =0 . 

T  D  ~  0  . 

QW=0  • 

BPOS=  1.0 
BMEG=-1 .0 
BD  =  BRB(  1,2) 

NK=MOD ( NORD ( I ) , 2 ) + 1 
CO=COS( P ( I ) ) 

C  F  =  C  0  S  (  R  (  I  )  ) 

S  D  =  S  I  N  (  P  (  I  )  ) 

SE=SIN( R( I ) ) 

TP  1  =  6,233185 
XK=TPI*RIO( I ) 

TF  =  EXP  (  -BR8  (1,1  )  *R  HOSQ )  *GP(  I  +  NU  )  *  I  C  ENT*NSTEP  (  I  ) 

IE (NAD. EO. 2)  TF=TF#2 • 

TF=TF*F( JF) 

DO  14  IWA=l,2 
DO  10  fW  =  l  ,6 


. 
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DEL  MCI  IWA, IW)=0. 

DELMSt I W  A , IW)=0.0 
IF ( IW.GT.3)  GO  TO  10 
OF  LC (  IWA,IW)=0. 

DELS! I W A  , IW)=0. 

10  CONTINUE 

14  CONTINUE 
CODE^O . 

on  ? 3  J  =  lfMU 
I  CC!DE  =  0 

CALL  SYM 

PHI=TP I*(RH*XC( I ) +RK*YC( I )+RL*ZC ( I } +T ) 

CP=COS( °H I )  *TF 
SP=$IN(PHI )*TF 
XK  l=RH*BAl+PK*8A2+RL*BA3 
XK2=RL*BA6 
X  K3  =  R  K*  B  A  4+ RL  *  8  A5 
C C  =  X K*  {  X K  1  * C E - X K 2*SE) 

D0=XK*(XKl*S0*SE+XK3*CD+XK2*Cr*SD) 

AP-=DS  OR  T  (  CC*CC+DD*OD) 

XNU-QARCOS  (CC/AP) 

I F ( DO  *L T . 0 • 0 ) XNU=-XNU 
ANGLE=NORD ( I ) * { U ( I J-XNU) 

I F { MODE . FQ . 1)  GO  TO  2 

I F { KSFL ( L+7 ) , EO.O . AND.KSELI L+3 ) . FQ. 0)  GO  TO  2 
DELCE=-XK* ( XK1*SE+XK2*CE) 

DELDE=XK*( XK1*SD*CE-XK2*S0*SE) 

0  F  L  0  D  =  X K * ( C D*S E*XK 1 -SD*XK3+CE* C  D* XK 2 ) 

AP2=AP*AP 

DEL APD= DD*DELOP/AP 
D  FLAP  F  = (CC*DELCE+DD*DELDE )/AP 
TE  =  DEL A  PE  ZAP 
TD=DELA  PD/AP 

******  THIS  NEXT  SECTION  WILL  CAUSE  DIVIDE  BY  ZERO  ERRORS 
FOR  CERTAIN  REFLECTIONS  IF  D  DR  E  ARE  0,90, ISO, OR  270DEGS 
BY  PASS  IS  TO  OMIT  THESE  REFLECTIONS  FROM  LEAST  SQUARES 

CW=AP2*DSQRT ( 1 .-CC*CC/AP2) 

IF  (ABS (CW) .LT. 0.00001)  GOTO  1 
DXNU=1. /CW 

DXNUE=-DXNU*{ AP*DELCE-CC*DFLAPE) 

DXNIJD  =  DXNU*  C  C  *  D  E  L  A  P  D 
GO  TO  2 

1  I  CODE= 1 

2  IF ( NK.EQ.2)  GO  TO  3 

*****  CALC 1 N  OF  REAL  TERMS  ***** 

15  MPRIME=NK*NOP.D(  I  )/2 
NA  A=  1 

GO  TO  21 
22  F  =  SNGL ( AP ) 

CALL  BESSEL  (BJ,0W,F  , BNEG ) 


. 


- 
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s  u  m ( i ) - n  b  r  b ( r )*bj 

SUM ( 2 )  =  BR 
SUM [ 3 ) = BOD 
SUM {  4  )  =  BE  A 
S  U  M  {  5  )  =  0  •  ' 

SUM(6)=BJ 
GO  TO  5 
C 

C  *****  C AL C 1 N  OF  IMAG  TERMS  ***** 

C 

3  NPR I HE~ { NORD  ( T )-l) /2 
DO  4  IA=1,6 

SUM(IA)=0. 

4  CONTINUE 
NA  A-2 

21  DO  24  IA=l,5 

LM( I  A )  =  K  S  E  L  ( IA  +  L  +  4) 

24  CONTINUE 
L  M  (  6 )  =  l 

IF(NAA.EQ.l)  GO  TO  22 

5  PS  I GN=  1  • 

IF  ( ICODE.FO.O)  GO  TO  16 
SUM(3)=0. 

SUM(4 )=0. 

L  M  (  3  )  =  0 
L  M  (  4  )  =  0 

16  DO  6  IA=1,NPRIME 

PSIGN=PSIGN*(-1) 

6  CONTINUE 
SIGNERS IGN 

IF ( NK . F  Q. 2 )  PSIGM  =  -1* 

NL=1 

IF ( NA A . ME .NK) NL=2 
C 

c  *****  CALC 1 N  OF  M  AND  DERIVATIVES  ***** 

C 

DO  11  JP=NL,20fNK 
J  PA- J  P 

PA=F LOAT { J  P) 

PNA=PA*NORD( I ) 

CU=COS( PA-ANGLE) 

SU=S IN (PA* ANGLE) 

TC  OM  =2 . *  SI GN*B I A ( I , J  P ) / B 1 0  (I) 

F  =  SNGL ( AP) 

CALL  BESSEL  (BJ,PNA,E  ,RNE3) 

IF {MODE. EO.  1)  GO  TO  34 

I  F ( L  M { 1 J . E 0 . 0 )  GO  TO  30 
C 

c  ******  CALC  *  N  OF  DM/D (3D)  ***** 

C 

3DT  E  R  M-  2  .  *  S  I  GN  *  C  U  *  DO  J  B  (  T  ,JP  )  *BJ 
SUM( i )  =  SUM ( 1 ) +BDTFRM 

I  F  (  APS  (  BDTFR/1/SUM  (  1)  )  .LE.CLIMIT)  L  M  (  1 )  =  0 
30  I F ( L M ( 2 )  .EQ.O)  GO  TO  31 


. 

■ 
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*****  CALC  *  N  OF  DM /OR  ***** 

R  T  E  R  M=  TC  CM  *  B  R  *  C  U 
SUM ( 2 )  =  S  UM ( 2 ) +  R  T  ERM 

I F ( A3S ( R TFRM/SUM ( 2 ) )  , LF . CL  I M  I  T )  LM ( 2 )  =  0 

31  r  F ( L M ( 3 ) . F Q . 0  )  GO  TO  32 

*****  CALC’N  OF  DEL M D  ***** 

D  M  D  -  S  U  *  P  NA  *  D  X  N  U  D 

DT  FR M=  TCOM * ( 3D0*CU  +  B  J  * DMO ) 

SUM( 3)=SUM(3 ) +OTERM 

IF { A  B  S ( D  T F  R  M / S  U  M ( 3 ) )  . LE.CLIMIT)  LM(3)=0 

32  IF(LM(4)  • EQ . 0 )  GO  TO  33 

*****  CALC 1 N  OF  OM/OE  ***** 

D  M  E  -  S  Li  *  P  N  A  *  D  X  N  U  F 
ETERM=TCOM* ( BE A*CU+B J*DME ) 

SUM ( 4 )  =  SUM ( 4 ) +ETFRM 

IF( ADS (ETFPM/SUM(4) ) , LE.CLI MIT )  LM{4)=0 

33  IF(LM( 5)  .  FQ . 0 )  GO  TO  34 

*****  CALC'N  OF  DM/DG  ***** 

GT  F R  M=  - P  M A  *T  C  OM  * S  U *R  J 
SUM!  5)=SUM(5)  +-GTFRM 

I F ( ABS ( G TERM/ SUM ( 5 ) )  . LE. CLI M I T ) L  M ( 5 ) =0 

34  IF(LM(6) . EQ . 0 )  GO  TO  35 

*****  CALC’N  OF  M  ***** 

MTERM=TCOM*CU*RJ 
SUM  {  6 )  =  SUM  (  6  )  4-MTERM 

I F ( A  BS ( M  TERM/SUM  16)). LE.CLIMIT)  LM ( 6 ) = 0 

35  L5UM-LM ( 1 ) +L M ( 2 ) +LM ( 3 ) +L M { 4 ) +LM ( 5)+LM(6) 

I F ( MODE  *  EQ .  1 )  LSUM=LM(6) 

Ic (LSUM.EQ.O)  GO  TO  12  . 

S I GN  =  S I G  N*  PS  I GN 

11  CONTINUE 

WRITE (6 ,201  )R!I,RK,RL 

12  IF  { MODE. EQ ♦ 1 )  GO  TO  7 

SUM! 1 )= ( SUM ( 1 )— SUM ( 6)*DBIB{ I )) /BIO{ I ) 

DO  34  I A=l , 5 

DELMC { N  A A , I  A) =  SUM(  I  A } *C P  +  OFl MS ( MAA,  T A) 

DFLMS ( NAA , IA)=SUM( I A)*SP+DELMS (NAA, IA) 

36  CONTINUE 

IF ( (KSELIL+l )+KSEL ( L+2) +KSEI I L+3) ) . EQ.O)  GO  TO  7 

OFLCIMAA,  1 )=OFLC (NAA,  1  )  +CP*SLM  {  6)*RH 

DE  LC ( NA  A , 2 ) =OELC ( NAA , 2 ) *C P*SUM ( 6) *RK 

OELC( NAA, 3 ) =OELC (NAA, 3) »-CP* SUM (6)*RL 

DELS (NAA, 1)=DEL S ( NAA ,  1)+SP*SUM ( 6 ) *R  H 

DFLS( NA A, 2) =0 FL S ( N A A , 2 ) +S P* SUM ( 6) *RK 


' 

' 


DH  L  S  (  N  A  A ,  3  )  =n  FL  S  ( ! '  A A , 3 ) +S P* SUM l  6 )  *R  L 
7  HE  I.  MC  (  NAA  »  ft  )  =  SUM  (  6  )  *C  P+  DEL  MC  {  M  A  A  ,  6  ) 
nELMS(NAA,6)  =  Sl)M(  6)  *SP  +  DELMS(  NAA,6) 

IF (NAA. EQ.2 )  GO  TO  15 
23  CONTINUE 
3  AT  =  DELMC<  1 , 6 ) -DELMS ( 2 , 6) 

AC AL=ACAL+AT 
I  F ( MAD . FQ . 2  )  GO  TO  40 
3  T  -  D  E  L  M  C  (  2  ,  6)  +0  F  L  M  S  (  t  ,6) 

BC AL=BCAL+BT 

40  IF ( MODE .EQ.  I  )  GO  TO  9 

*****  CALC'N  OF  D F L F (POPULATION  FACTOR)  ***** 

IFlKSEL(l) .EQ.O)  GO  TO  58 
DA ( K ) = AT/GP ( I+NU) 

IF  C INV1  .FQ.  I )  GO  TO  59 
DB(K)=3T/GP( I+NU) 

59  K=K+1 

53  L=L+ 1 

*****  DERIVATIVES  W.R.T.  CENTRE  OF  RING  ***** 

DO  56  I A=  1 , 3 

IF ( K S E L ( L ) .EQ.O)  GO  TO  55 
D  A ( K ) =— T  P I  * (DELS l 1 »I AJ+PELC (2, I  A)  ) 

IF ( IMV1 .EQ. 1 )  GO  TO  57 

DB (K )=TPI* (-DELS (2,1  A) +DELC ( 1,  I  A  )  ) 

57  K  =  K  +  1 

55  L=L+1 

56  CONTINUE 

IF (KSEL (L) , EQ.O)  GO  TO  53 

*****  C A L C 1 N  OF  DEL F (TEMP)  ***** 

DA ( K ) =-RHOSQ* ( DELMC ( l, 6) -DEL MS (2,6)) 

I F ( INV1  .EQ. 1 )  GO  TO  54 
DB(K)=-RHOSQ*( DELMC (2,6) +DELMS (2,6) ) 

54  K=K+1 
53  L=L+ 1 

*****  CALC  *  N  OF  DELFD  TO  DFLF3  ***** 

DO  50  I  A=  1 , 5 

IF  (KSEL.  (  L)  .EQ.O)  GO  TO  51 
DA (K  )  =  D E L M C ( 1, T A ) -DELMS ( 2 , I  A) 

IF ( I NV1. ED. 1 )  GO  TO  52 
DB(K)  =  DELMC(2, I  A ) +DELMSI 1 , I  A ) 

52  K=K+ 1 

51  L  =  L+ 1 

50  CONTINUE 
9  RETURN 

201  FORMAT ( *  REFLECTION  » ,3F3. O,1  HAS  NOT  CONVERGED*) 
END 


■ 


' 

' 


■ 
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C 


C 


SUBROUTINE  BESSEL  (BSIJM,  PP,  X7.  ,  BSIGN  ) 

CO MMON  T  I  TL  E  (  1 8  )  ,  M.ODF  ,  I  MV , NAD,  I S A  N , NP  J , NEW ( 5 )  ,  NF SQ , NE  OUR , 

1  NO  OR  ,  KAR  D  ,  I  0 ,  L  V  ,  NLV  ,  NS  ,  MU  ,  MU  ,  M  LSC  ,  K  R I T  ,  P.P  AR  ,  N  AT  E,  IC  EN  T 

2  ,  N  R  I G  ,  MDNR  ,  M! !  I  N 

COMMON  A,B,  C, AL, BE ,GA , AQ, BQ,C3  , AL Q , BEQ , GAQ , NSL , SUMDL 
COM  MON  NP AR , N VAR , NC  NT , K  SEL ( 921 ) 

COMMON  S  (  50  )  ,  BG  ,  GP  {  70  )  , GQ ( 60 ) , GR ( 60 ) , ME ( 60) , MG ( 60 ) , M I { 60 ) , 
2X{ 140) ,Y{ 14  0) i Z { 1 40  )  ,R1  l(  140),  822(60) ,B 33(60)  ,B 12(60) » 
2813(60)  ,028(60)  , N A ( 7 0 )  , NAME (140)  , X C (  i 0 )  ,YC(  10)  ,ZC(10)  , 

3P(  10) ,R (10)  ,U( 10) ,8R8( 10,8) ,N>TEP(  10) , NGRD (  10)  , R l 0 (  10) 

COMMON  MLSVJ,  AW  ,  BW  ,  C  W ,  DW  ,  EW  ,  E  W ,  GW,HW  ,  S  IZE1  ,S I ZE2 
COMMON  SAER (  10,  10 )  , NSAER (  10 , ID )  , SAE S ( 10 )  , NS  AES ( 10 ) , SAFE ( 10) , 
lMSAEF(lO) 


COMMON  NTR AN ( 4 , l 2 , 1 2 ) , N VTN ( 4, 3 6 ) 

COMMON  NE , E M ( 3 2 , 20 )  ,F(20)  ,EP(?0) ,FPP(20)  ,RHO,RHOSQ 
COMMON  FOBS , FCAL , S  FC , AOBS , ACAL , BOBS , BC AL , DEL , S  IGMA , EXT l ,E  XT2 
1  t  M  S 1  ,  MR  E  J  ,  Q  ,  MH  ,  MK  ,  »ML,  TH  ,  T  K  ,  TL  ,  D  A  (  180), DB(  130)  ,  DF  (  1 3  0  ) 

COMMON  C  H  G ( 0  2  1  )  ,STD(921 ) ,EPS(130)  , D I  AG ( 130)  ,V(  180) 

COMMON  All, A] 2 , A  1 3 ,A22, A23,A33 ,BAl ,BA2,3A3,BA4,RA5,BA6 
COMMON  SN( 50)  ,BON,GPN( 70)  ,CQN( 60)  ,GRN(60)  , XM ( 1 40 )  , Y N ( l 40 )  , 
17.N14C)  ,  BUM(140),B22N(60)  ,P33vj(  60),B12N(60)  ,B13N(60)  , 

2B23N( 60) ,XCN( 1 0) ,YCN( 10 ) , ZCNd 0 ) , PN ( 10 ) , RN ( 10 ) ,UN ( l 0) , 

3  B  R  B ( 10,8) , R IDN ( 10 ) 

COMMON/ SPSY / I , J , N , RH, RK , R  L , T, R HH , RKK , RLL , RHK , RHL , RKL , NR 
COMMON  /DDC/SPK,  SPL  ,SPP  ,SPP.,SPJ  ,CPH  ,  CPK  ,  CPL  ,  CPP  ,  CPR, CPU, WL  , 

1  SACAL,  W3, SBCAL,K, L,TPT , JF ,LC  ,SPH,INV1 
COMMON/ BE  S  S / OBJ  8 ( 10 , 20 )  , LM ( 6 ) ,  BR , BE  A, TE , TO, BDD , 3D , B 1 A ( 10 , 20 
1, DBIR (10), 310(10), J PA 
COMMON/ RIG  I O/CL IM IT 
DIMENSION  W L (20) ,WR(20) 

EQUIVAI.  EMCE  ( B  I ,  B 1 1 )  ,  (  B I  N ,  B1  IN)  ,  (  AR,  SN) 

01  MENS  I  ON  B I ( 14 0 ) , B I N ( 1 40 ) 

DIMENSION  A R (20100) 

DIMENSION  LB (5) 


BTERM=l . 

RRTER  M=0. 

FRTERM=0. 

DRTEP.M=0. 

BJTERM=0. 

L  B  {  5  )  =  l 

IF ( PP.EQ.O. )  GO  TO  17 
FACP=1. 

J  P -  I E I  X (PP) 

no  l  m= l » jp 

XM= FLOAT { M ) 

BTERM=  BT  EP.M*XZ/ ( 2 .*XM) 

1  CONTINUE 

17  IF  I BS IGN.LE.0 . )  GO  TO  2  2 

*****  CALC ’ N  OF  I  TERMS  ***** 
LB (  1 ) =1 
DO  13  M  =  7  »  4 
LB  f M )=0 
13  CONTINUE 


■ 


■ 
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C 


?  2 


14 


5 


in ( PP.NE.O. )BJTFRM=PP*BTFRM/B:> 


D B J 3 { I t  JPA) =BJTERM 
GO  TO  5 

*****  CALC  *  N  nr  H  TERMS 
LB ( 1 ) =0 


DO  14  M  =  2 1 4 
L  B  (  M )  =  L  M  ( M  ) 


CONTINUE 

IF (PP.EQ.0. )  GO  TO  5 
BR  T  ER  M=  BT  ER  M*  P  P  /  R  I  D  (I) 
ERTEPM=PP*BTERM*TF 
ORTERM=PP*BTERM*TD 
BSUM= BT  ERM 


BR=BRTE RM 
BE  A=  ER TERM 
BDD=DRT  ERM 

I F ( 3S i GN. GT . 0 }  OB J B ( I f  J  PA )  =  R JT  ERM 
00  2  M=  1 , 30 

BTERM=(BSIG N  *BTF  R M *XZ*XZ)/(4. * M * ( M  +  PP )  ) 
IF (MOOE.EQ. 1 )  GO  TO  11 
I  F  (  LB  {  1  ) . FQ . 0 }  GO  TO  6 


*5 


nI/  nV  »•/ 
V  'i'  'r 


CALC 1 N  OF  DM/ DSD  TERM 


>)c  j{e  jjc 


B J  TERM= ( 2 • *M  +  PP ) *BTFRM/BD 

OBJS{ I , JPA) =DBJB( I T  JPA) +BJTERM 

IF ( ABS( BJTFRM/DBJB ( I , JPA)  ) . LE. CL  I  MI T) LB( 1 ) =0 

6  IF  { LB  C 2) •FQ.O)  GO  TO  7 

*****  CALC 1 N  OF  DM/DR  TERM  ***** 

B  R  TER  M=  BT E R  M*  (  2  •  * M  +  P P  )  /  R I  D  (  I  ) 

BP  =  BR +BRTERM 

I F { ABS ( BRTERM/BR ) . LF . CL  I M  I  T )  L  B ( 2  )  =  0 

7  IF  ( LB ( 3) .EQ.Q)  GO  TO  8 

*****  C  ALC  *  N  OF  DM/DD  ***** 

D  R  T  ERH=  B  T  ER M * ( 2 . * M  +  PP ) *TD 
BDD=BDD+QRT  ERM 

IF  t ABS( DRTFRM/BDO) .LE.CLI MIT)  LB(3)=0 
3  IF  ( L  B ( 4 ) .FQ.O)  GO  TO  11 

*****  CALC 1 N  OF  OM/OF  TERM  ***** 

ERTERM= BTER M* (2 .*M+PP ) *TE 
BE  A=  BE A  +  CRT ERM 

IF ( ABS( ERTERM/BEA) . LE .CL  I  MIT)  L  B ( 4 ) =  0 


. 


■ 
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11  IF  l  LB  t  5)  .  f: Q  «  0 )  GO  TO  1  ?. 

*****  CALC ' N  OF  BFSSEL  SUMMATION  ***** 
BSUM---BSUM  +  BTFRM 

I F (  ABS ( BTERM/BSUM ) .  LE.CLm  T)  LB  15) =0 

1 2  L  SUM=LB (1)+LB(2)+LR(3)+LB(4)+LB(5) 

I F ( MODE  *FQ . I )  L  SUM=LB { 5 ) 

IF(LSIJM.EQ.O)  GO  TO  3 

2  CONTINUE 

3  RF TURN 
FNO 


X 
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PROGRAMMES  USED 
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Author 

Title 

Description 

D.P.  Shoemaker 

MIXG2 

calculates  Picker  diff¬ 
ractometer  settings 

modified  by  M.  Elder 
to  sort  data 

M.J.  Bennett 

PMMO 

calculates  intensities, 
makes  Lp  corrections  for 
Picker  data 

M.  Elder 

Datap 

as  for  PMMO,  except  for 
PAILRED  data 

M.  Elder  and 

K.  Simpson 

D-ref ine 

refines  axial  lengths 

A.  Zalkin 

FORDAP 

Fourier  summation  for 
Patterson  or  Fourier  maps 

M.J.  Bennett 

MMMR 

calculates  starting  para¬ 
meters  for  rigid  bodies 
and  hindered  rotors 

W.C.  Hamilton 

G0N09 

absorption  corrections  - 
Picker  data 

D.P.  Shoemaker 

GNABS 

absorption  corrections  - 
PAILRED  data 

C.T.  Prewitt 

SFLS5HR 

structure  factor  calcu¬ 
lation  and  least  squares 
refinement  of  parameters 

modified  by  B.M.  Foxman, 
M.J.  Bennett  and  W.L. 
Brooks  for  rigid  body  and 
hindered  rotor  routines 

G.J.  Williams 

CROMERS 

calculates  form  factor 
curves  using  Cromer 
coefficients 

J . S .  Wood 

MGEOM 

calculates  bond  lengths, 
angles  and  best  planes 

. 

" 


' 
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Programmes  Used 
Author 

W.  Busing  and 
H . A .  Levy 

C.  Johnson 

M.  Elder  and 
W.L.  Brooks 


(Contd. ) 

Title 

0RFFE2 


ORTEP 

SFLIST 


Description 


calculates  bond  lengths, 
angles  and  associated 
errors 

modified  by  B.  Penfold 
for  IBM  360,  and  W.L. 
Brooks  and  M.  Elder  for 
hindered  rotors  and  rigid 
bodies 

writes  plot  command  tape 
for  Calcomp  plotter 

lists  structure  factors 
for  publication 


. 


Appendix  C 


CONVENTIONAL  CRYSTALLOGRAPHIC  SYMBOLS 
AS  DEFINED  IN  VOL.  I,  PAGE  xi  ,  OF  THE 
INTERNATIONAL  TABLES  FOR  X-RAY  CRYSTALLO¬ 


GRAPHY 
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h  ,k ,  i 

indices  of  the  reflection  from  a  set  of 

parallel  planes;  co-ordinates  of  a 

reciprocal  lattice  point 

F 

hk£ 

structure  factor  for  the  unit  cell, 

corresponding  to  the  Bragg  reflection 

hk£ 

a,b,c. 

lengths  of  unit  cell  edges 

a  ,3  ,Y 

interaxial  angles 

a* ,b* ,c* 

lengths  of  reciprocal  lattice  unit  cell 

edges 

a  * , 3  *  *  Y* 

interaxial  angles  in  reciprocal  space 

x .  ,y .  ,  z . 

i  ’ 1  l '  l 

fractional  co-ordinates  of  an  atom  i 

(co-ordinates  of  atom  i  in  units  of 

a,b,c) 

p 

x,y ,  z 

electron  density  at  the  point  x,y,z 

y 

linear  absorption  coefficient 

-2 

u 

mean  square  amplitude  of  atomic  vibration 

B 

Debye  isotropic  temperature  parameter 

_  2-2 

B  =  8  TT  u 

311'322'333' 

anisotropic  temperature  parameters  used 

ei2'ei3'623 

to  describe  ellipsoidal  electron  distri¬ 
bution  of  the  anisotropically  vibrating 

atom;  the  temperature  factor  expression 
is  then:  exp  [- ( $ ^h^  +  322^2  +  ^33^^ 

+  2  3 12hk  +  2313h£  +  2323k&)] 
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F 

o 


F 

c 


observed  structure  factor 


calculated  structure  factor 


